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ABSTRACT 


In  work  related  to  radiation  shielding,  the  use  of  radioisotope  techniques, 
and  activation  analysis,  an  experimenter  rr  .  t  often  analyze  counting  data 
where  counts  are  caused  by  the  natural  background  and  by  the  decay  of  more 
than  one  radioisotope.  In  this  report  a  procedure  is  developed  for  estimating 
the  strength  of  each  isotope  at  different  times  from  several  decaying  radio¬ 
active  samples  of  a  single  multiple-isotope  source.  I  n  addition,  the  procedure 
provides  a  method  for  placing  confidence  limits  on  the  strengths  and  a  method 
for  partitioning  the  imprecision  of  estimating  the  strengths  into  thme  principal 
causes:  Poisson  variation,  sampling  error,  and  residual  error  (called  mechanical 
errcr). 

An  operational  FORTRAN  ll-D  computer  program,  SAND,  implements 
the  procedure.  The  procedure  and  program  were  tested  by  using  fictitious 
data  with  known  properties  as  inputs.  The  results  of  the  simulation  were  in 
reasonable  agreement  with  the  theoretical  values. 
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INTRODUCTION 

In  various  nuclear  physics  laboratories  there  exist  dilute  solutions  of 
mixtures  of  a  few  radioactive  isotopes  for  which  an  estimate  is  desired  of  the 
expected  rate  of  radioactive  ducay*  in  counts  per  minute  per  given  aliquot 
quantity  of  the  solution  at  some  given  historic  instant  of  time.**  Usually  it 
is  desired  to  break  this  total  estimate  into  the  parts  associated  with  each 
isotope. 

To  accomplish  this,  it  is  customary  to  take  several  aliquot  samples  and 
observe  counts  from  each  of  these  during  measured  short-time  periods  at 
several  different  ages  (lengths  of  time  from  the  desired  historic  instant  t„  *** 
to  the  beginning  of  the  measured  period).  Because  many  materials  in  nature, 
and  particularly  in  a  nuclear  physics  laboratory,  are  experiencing  some  radio¬ 
activity,  the  observing  device  (Geiger  counter)  will  count  slowly,  even  v.ith 
no  sample,  thus,  it  is  necessary  to  record  counts  of  this  background  during 
relatively  long,  measured  time  periods  near  the  time  when  the  counter  is  to 
be  used  for  sample  counting.  The  purpose  of  this  background  count  is  to 
correct  the  observed  counts  to  obtain  apparent  sample  counts. 

Given  all  the  foregoing  data,  the  following  questions  arise: 

1.  At  the  desired  historic  instant,  designated  ,  what  is  the  estimate  of  decay 
rate? 

2.  At  t^  what  is  the  breakdown  of  the  total  decay  rate  into  the  various  isotope 
sources? 

The  laboratory  data  used  to  make  the  foregoing  estimates  will  be  subject 
to  the  following  three  sources  of  imprecision: 

a.  Sampling  error;  that  is,  inability  to  place  exactly  the  same  amount  of 
material  in  exactly  the  same  position  on  the  watch  glass  each  time  in  each 
finally  prepared  sample  ready  for  counting  and  recounting. 

*  The  terms  "rate  of  decay”  and  "strength  at  some  instant  of  time"  are  used  inter¬ 
changeably  throughout  this  report. 

**  U.S.  Naval  Civil  Engineering  Laboratory.  Technical  Note  N-559:  Interim  calibration 
and  procedures  for  GM  detector,  by  L.  B.  Gardner  and  B.  Burdick,  Port  Hueneme, 

Calif.,  Mar.  1964. 

***  The  reader  is  referred  to  the  Glossary  of  Symbols  on  a  holdout  page  at  the  back  of 
this  report. 
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b.  Poisson  distributed  variation;  that  is,  the  fact  of  nature  that  a 
multiplicity  of  samples  of  exactly  the  same  size  and  strength  (same  expected 
counts  per  minute)  will  show  variation  among  themselves  in  the  actual  number 
of  radioactive  events  in  a  given  minute.  The  nature  of  this  Poisson  fluctuation 
is  known. 

c.  All  other  uncorrected  sources  of  imprecision  lumped  under  the 
arbitrary  heading  "mechanical  variation."  These  sources  may  include  slight 
positioning  variation,  extraneous  counting  errors  other  than  from  background, 
temperature  variation  of  mechanical  setup,  slight  gain  or  loss  of  material  in 
"dried"  sample  between  successive  observations,  and  inexactness  of  dead¬ 
time,  r,  correction  used. 

Because  of  a,  b,  and  c  the  estimates  1  and  2  will  not  be  exact.  This 
raises  the  further  questions: 

3.  How  imprecise  is  the  estimate  of  total  decay  rate  at  t'  ? 

4.  How  imprecise  is  the  estimate  of  each  isotope  decay  rate  at  ? 

An  incidental  result  of  attempting  to  supply  an  answer  to  3  will  be  an 
approximate  answer  to 

5.  How  much  of  the  imprecision  in  estimating  decay  rate  at  is  due  to: 

(a)  sampling  error,  (b)  Poisson  error,  and  (c)  mechanical  error? 

It  is  the  purpose  of  this  report  to  develop  a  computational  procedure 
for  approximating  answers  to  1 , 2,  3.  4,  and  5.  In  short,  these  five  questions 
constitute  the  problem  to  be  solved  in  this  report. 

Briefly,  the  computational  procedure  will  consist  mainiy  of  a  sequence 
of  parallel-weighted  linear-least-squares  curve  fits.  On  each  iteration  each 
datum  is  weighted  inversely  proportional  to  its  current  estimated  variance. 

The  repetitions  are  terminated  when  the  last  two  current  estimates  of  the 
mechanical  error  standard  deviation  differ  by  less  than  some  arbitrarily  small 
value.  The  sampling  error  standard  deviation  is  estimated  by  an  additional 
iterative  procedure  using  the  foregoing  results. 

It  is  assumed  that  the  number  of  isotopes  in  the  solution  and  their  decay 
constants  are  known;  that  is,  the  procedure  does  not  estimate  the  number  of 
isotopes  in  the  solution,  nor  does  it  estimate  the  decay  constants  since  better 
methods  exist  for  this  purpose. 


DESCRIPTION  OF  MATHEMATICAL  CONCEPTS 
Ganeral 

Several  aliquot  samples  (total  of  I  samples)  are  taken  from  a  liquid 
solution  of  radioactive  material,  placed  on  I  watch  glasses,  and  allowed  to 
evaporate.  Then  each  sample  is  observed  several  times  by  means  of  a  detector 
(J,  times  for  the  ith  sample,  i  =  1,2,  I)  in  a  fixed  counting  system  with 
fixed  geometry.  The  number  of  counts  recorded  during  a  period  of  duration 
Tjj  minutes  will  be  designated  C|j  for  the  i^  sample  and  its  jt*1  observation 
(i  =  1,2, I;  j  =  1,2, Jj).  The  age  of  the  sample  in  minutes,  .eckoned 
from  t' ,  an  arbitrary  historic  instant  of  interest  when  t  =  0,  to  the  beginning 
of  the  period  of  duration  T jj,  will  be  designated  t^. 

The  assumed  known  number  of  radioactive  isotopes  in  the  source 
solution  is  L.  The  known  decay  constant  in  reciprocal  minutes  from  the  8^ 
isotope  is  designated  Xg(8  =  1,2,...,  L). 

Because  the  detector  records  counts  from  background  radioactivity  in 
addition  to  those  striciiy  from  the  sample  under  observation,  an  effort  is 
made  to  correct  for  this  by  observing  background  with  the  same  instrument 
at  a  time  near  t,j.  The  length  of  this  observation  and  the  number  of  counts 
observed  are,  respectively,  T,*  and  C,*.  In  practice  "near"  may  mean  merely 
in  the  same  working  day.  Thus,  a  single  background  reading  may  be  used  to 
correct  several  observations  and  give  rise  to  T^  and  C|)  combinations  which 
are  numerically  identical  for  various  i  and  j  combinations. 

In  the  sampling  process  an  effort  is  made  to  place  the  same  amount  of 
material  on  each  watch  glass  in  the  same  relative  position.  There  is  a  slight 
variation  in  this  process,  and  the  effective  quantity  will  be  assumed  to  be  a 
normally  distributed  random  variable  with  mean  Q  and  standard  deviation 
Qot  so  that  os  may  be  construed  as  a  proportional  standard  deviation  of  error. 
For  example,  at  =  0.01  would  indicate  a  sampling  error  standard  deviation  of 
1%  of  the  expected  effective  sample  size. 

A  period  of  time  elapses  between  successive  observations  of  the  same 
sample,  a  period  in  which  laboratory  conditions  change.  An  effort  is  made 
to  make  all  observations  with  exactly  the  same  mechanical  and  electrical 
arrangement.  Nevertheless,  there  may  be  a  slight  variation  in  this  process 
which  affects  the  actual  Cjj  recorded.  If  one  could  imagine  a  fictitious  sample 
(sample  number  h)  which  would  experience  exactly  the  same  number  of 
radioactive  events  every  minute,  regardiess  of  age,  then  the  successive  obser¬ 
vations  ^hl  >  Ch2  -  Ch  ...  with  fixed  Thj  and  no  background  noise  will  be 
assumed  to  be  items  of  a  normally  distributed  random  variable  with  mean  M 
and  standard  deviation  M  om .  Now  am  may  be  construed  as  a  proportional 
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error  due  to  all  causes  (called  mechanical)  other  than  sampling  error  and 
errors  that  result  from  the  Poisson  phenomenon  (assumed  independent  of 
mechanical  errors). 

If  one  now  imagines  an  infinite  number  of  samples  of  the  same  solution 
of  radioactive  material  and  an  infinite  number  of  observations  of  each,  all 
between  t£  and  t'0  +  T,  for  fixed  T,  then  the  background-corrected  counts 
per  minute  would  constitute  a  two-dimensional  array  of  data  entries,  designated 
here  as  t}|(,  infinite  in  length  and  infinite  in  width.  The  row  (sample)  mean  of 
such  a  fictitious  array  of  data  will  be  designated  m.|(0,  T)  counts  per  minute. 
The  portion  of  this  due  to  the  £^  isotope  will  be  designated  Mg ,  (0,  T),  where 

M.|(0,T)  =  £  M8i(0,T) 

£-1 


If  the  average  of  the  infinity  of  the  m  t  (0,  T)'s  is  designated  as  (0,  T), 
then  question  1  of  the  Introduction  asks  for  m° ,  an  estimate  of  m°  ,  where 
zero  bias  error  is  assumed  and 

M°  =  lim  m  (0,  T) 

T-.0 


Similarly,  question  2  asks  for  the  L  Mg  's,  which  are  estimates  of  the  jig  's, 
where 


Mg  =  lim  Mg  (0,T) 

T-.0 


The  mean  of  the  i^1  sample  infinite  row  "at  t  =  0"  was  designated  as 
M.j  (0,  T),  where  it  is  supposed  that 

i  J' 

M.i  (0,  T)  =  lim  f  £  rji, 

V°°  1  j-i 

and  that  M  (0,  T)  =  lim  4-  V  m  i  (0,  T) 

i-~  i-i 

Ely  an  observation  "at  t  =  0"  is  meant  an  observation  of  small  duration 
containing  t  =  0  such  that  the  expected  counting  rate  at  t  =  0  is  equal  to  the 
average  expected  counting  rate  during  tne  time  interval  T.  Recall  that  this  is 
merely  a  fiction  inserted  to  clarify  concepts. 
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Briefly,  the  problem  is  to  obtain  estimates  of  m°  ,  m°  .  . m°  ,  o, ,  and 

om ;  to  compute  confidence  intervals  expressing  the  imprecision  in  knowledge 
of  the  values  of  Mi.. M2..  -  M°..  and  M°.  .'  and  finally  to  separate  the  imprecision 
in  estimating  m°.  into  the  three  sources:  (a}  sampling  error,  (b)  Poisson  variation, 
and  (c)  mechanical  variation.  To  per  form  these  estimate  -  it  will  be  necessary  to 
start  the  computation  with  omo ,  Xg,  t^,  T,j,  Cjjr  tjj,  T*,,  and  C-j,  where  6=1, 

2, ...,  L;  j  =  1,2, ....  J| ;  and  i  =  1,2 . I.  It  should  be  carefully  noted  that  this 

whole  report  is  concerned  solely  with  imprecision  and  not  laboratory  mean 
bias  over  the  testing  period.  Any  such  bias  in  m°  will  be  carried  over  into  m"  , 
and  as  such  will  remain  unkown  in  this  computation. 

To  reach  a  "solution"  to  the  problem,  least-squares  estimates  of  m°i .  .... 

Mij  for  each  sample  will  be  found  by  using  weighted  linear  regression.  For  a 
minimum-variance  estimate,  the  weight  used  for  each  datum  should  be 
inversely  proportional  to  the  variance  in  knowledge  (ignorance)  of  the  "true" 
(biased  by  the  particular  sample  bias)  value  represented  by  the  datum.  Given 
a  particular  sample,  this  variance  for  each  datum  is  the  sum  of  those  caused  by 
Poisson  phenomenon  and  mechanical  variation.  (The  Iasi  two  statements  taken 
together  contain  a  r,ight  inexactness  stemming  from  a  small  interdependency 
between  those  elements  of  data  which  were  computed  using  a  common  back¬ 
ground  reading.) 

Ideally  the  weight  of  a  datum  should  be  based  upon  the  true  variance  for 
that  datum;  however,  only  an  estimate  of  the  variance  is  available.  The  curve¬ 
fitting  method  depends  upon  the  weights,  which  in  turn  depend  upon  the  fitted 
parameters.  Thus,  it  is  necessary  to  iterate  on  weights  and  fitted  parameters 
simultaneously  for  all  samples  to  convergence  of  the  estimated  mechanical 
error.  This  leads  to  a  stable  om  and  I  stable-weighted  least-squares  curve  fits. 

Preliminary  Theoretical  Discussion 

In  the  foregoing,  the  word  "expected"  has  been  used  in  its  technical  sense 
and  as  such  does  not  mean  "anticipated."  For  example,  consider  an  infinitely 
large  quantity  of  radioactive  solution  and  the  division  of  a  large  portion  of  this 
with  no  sampling  or  mechanical  error  into  a  large  number  I  of  exactly  equal¬ 
sized  subquantities.  Now  imagine  I  background-free  detector  setups  with 
absolutely  no  variation  between  the  setups.  At  the  timet,,',  initiate  simulta¬ 
neously  the  I  countings  (one  subquantity  for  each  counter)  and  terminate  them 
all  simultaneously  T  minutes  late1 .  This  gives  I  integer  counts  C,  (i  =  1 , 2, ...,  I). 
Let 


e  -  {  £  c. 

i-1 
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If  I  is  allowed  to  increase  indefinitely,  C  will  approach  a  stable  value  known 
as  the  expected  value  of  C.  Restated 

E(C)  =  lim  C 


Thus,  the  expected  value  of  a  random  variable  is  merely  the  average  of  the 
numerical  outcomes  of  an  infinite  number  of  trials.  The  infinite  number  of 
outcomes  is  the  population  associated  with  the  random  variable.  A  few  out¬ 
comes  are  a  sample  of  the  population.  In  general  the  individual  C/s  will  show 
considerable  variation  among  themselves.  This  variation  of  the  population  of 
counts  has  been  called  the  Poisson  phenomenon.  Indeed,  the  a  priori  proba¬ 
bility  that  one  particular  counter  will  count  exactly  h  radioactive  events  during 
the  period  t'  to  +  T  is  given  by  the  Poisson  formula 


P(C,  =  h) 


e-E(c>  [E(C)]h 
h! 


for  h  =  0,  1,  2, ... 


Now  if  the  variance  of  C  is  defined  to  be 

oc2  =  E  j  [C  -  E(C)]  2| 


it  will  be  found  that 


ac2  =  E(C) 

or  equivalently  that  the  standard  deviation  of  counts 

oc  =  +  yjEiC) 

which  is  a  very  interesting  property  of  a  Poisson  distributed  random  variable 
such  as  C.  Restated 


oc2  =  [Cj  -  E(C>] 2 

'■*"  i=i 

Because  counts  per  minute  for  no  detector  deadtime,  r  =  0,  is  N  =  C/T,  T 
merely  plays  the  role  of  a  scale  factor  with  the  following  consequences: 
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Of  course  in  a  real-world  laboratory  it  is  not  possible  to  perform  the 
idealized  experiment  of  this  section  because  there  may  be  only  one  detector. 
In  such  a  case  several  samples  may  be  observed  at  different  ages,  or  each 
sample  may  be  reobserved  at  different  ages.  For  a  fixed  period  length  of  T 
minutes  the  expected  counts  wili  decrease  or  decay  with  age.  Fortunately, 
the  physical  law  governing  this  decay  of  expectation  is  known.  If  at  time  t' 
the  instantaneous  expected  background-free  counts  per  minute  are 

£  -  i  4. 

C"1 

then  at  any  future  instant  t 

km  •  £  4  e'*8* 

8  =  1 

and  over  a  period  from  t  to  t  +  T  the  expected  counts,  assuming  r  =  0  (perfect 
detector),  are 


E(C) 


P.'.(t)dt 


+  E(C*) 


where  C*  is  the  actual  counts  from  background  during  T. 

Equivalently,  the  expected  counts  per  minute  solely  from  sample  are 
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Letting 


E(tj)  =  j 


H'„  (t)  dt 


-Xp  t 
e  * 


then 


X(t,T;Xg) 


e(t?)  =  Y.  4  x^V 

C-1 


For  the  jth  observation  (j  =  1,2 . J  j)  of  the  i1*1  sample 


E<*>ii)  =  E  #4  X(tjj, T jj ; Xg ) 

C-i 


and 


"  E  *4  X(tij-Tij;XC>  +  ei) 

C-i 


where  e(j  is  a  random  error  consisting  of  the  algebraic  sum  of  the  random 
Poisson  sample  error  (reckoned  from  expected  counts  per  minute),  the  random 
Poisson  background  error  (reckoned  from  expected  background  counts  per 
minute),  and  the  random  mechanical  error. 

Thus  the  estimates  of 


E<»„>  *  t  “l  *8,i 

8-1 


can  be  obtained  from  a  weighted  least-squares  fit  of  L  +  1  columns  of  data 
(each  Jj  long).  The  column  headings  are  r\ X1j(  X2j, Xu  and  the  corre¬ 
sponding  entries  are  tj..,  X(ti(,  T  ;  Xp)  =  Xp..,  where  j  =  1 , 2, ...,  J.  and  C  =  1, 

2 . L.  Each  of  the  Jf  lines  of  data  will  be  weighted  at  each  iteration  inversely 

proportional  to  the  last  iterated  estimate  of  E{ejj 2 ).  After  studying  the  next 
section,  the  careful  reader  will  note  that  the  foregoin^  ''entence  is  very  slightly 
inexact.  Because  a  single  background  reading  may  be  ,ed  for  more  than  one 
combination  of  i  and  j,  the  next  section  deliberately  modifies  the  method  of 
weighting  so  as  to  reckon  each  background  count  once  in  the  whole  data- 
weighting  computation  system. 
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Clearly,  one  will  never  be  able  to  perform  enough  countings  to  display 
the  whole  population  of  outcomes,  but  must  be  content  to  perform  several 
tests,  giving  merely  a  sample  of  several  numerical  valued  outcomes  from  the 
population.  The  average  of  the  sample  of  outcomes  is  called  a  sample  mean 
and  as  such  constitutes  an  estimate  of  the  population  mean.  For  example 

c  =  eTc) 

In  general  this  is  the  best  available  estimate  without  further  testing.  Nonethe¬ 
less,  it  may  be  a  poor  estimate.  It  is  customary  to  quantify  the  poorness  or 
goodness  of  the  estimate  by  computing  a  confidence  interval  at  some  expressed 
level  of  confidence  for  the  quantity  estimated.  A  0.95  interval  for  E(C)  is  a 
range  of  possible  values  with  specified  lower  and  upper  end  points  where  the 
interval  has  the  following  interpretation:  “The  true  value  of  E(C)  is  not 
known;  however  we  are  95%  sure  that  whatever  its  value  may  be  it  does  lie 
in  this  range."  Thus,  as  the  precision  of  the  estimate  increases,  the  confidence 
interval  decreases  in  length.  In  this  report  answers  to  questions  3  and  4  of  the 
Introduction  are  expressed  as  0.95  confidence  intervals. 


PROBLEM  SOLUTION 

The  broad  plan  for  solution  will  consist  of  four  main  steps: 

1 .  All  I  of  the  are  computed.  Each  of  these  computations  uses 
the  J;  observations  associated  with  the  particular  sample;  amo,  the  historic 
estimate  of  proportional  mechanical  error  standard  deviation  with  f  associated 
degrees  of  freedom,  and  a  least-square  fit  to  the  X  weighted  points  to  arrive 
at  the  estimates  of 


L 

o  _  o 

M.i  -  D  Hi 

where  it  is  assumed  that 

*ig.  =  lim  .  (0,  T)  and  /u°.  =  lim  p,(0,T) 

T-0  T-0 

Estimates  of  om  based  upon  the  sum  of  historic  and  current  degrees  of  freedom 
are  used  in  the  computation  for  weighting.  To  obtain  om  all  I  samples  are  used. 
Thus,  all  I  least-square  fits  are  performed  in  parallel  using  the  results  of  the 
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0 


iteration  to  obtain  all  the  jk  (I  times  L  of  these)  and  omk  to  be  used  in  the 
(k  +  1)tf*  iteration.  The  magnitude  of  og,  the  proportional  sampling  error 
standard  deviation,  in  no  way  influences  these  computations  because  each 
fit  is  for  a  fixed  single  sample. 

2.  The  I  final  ?°.  's;  of0,  the  historic  estimate  of  at  with  d  associated 

degrees  of  freedom;  and  S  *,  the  estimated  variance  of  n°. ,  are  used  to  com- 

M°|  • 1 

pute  a  new  os  by  iteration. 

3.  The  statistics  from  the  foregoing  final  iterations  are  used  for 
computing  confidence  intervals,  and  for  partitioning  total  uncertainty  in 
?°  into  its  three  basic  components. 

4.  The  final  om  and  o5,  based  on  both  current  and  past  data, 
together  with  their  respective  associated  degrees  of  freedom,  are  retained  as 
inputs  to  future  data  analysis. 

Estimation  of  ^ .  and  om 

The  linear  mathematical  model  is 

V  =  E  Hi  xCij 

fi-i 

where  the  transformed  independent  variables  are 


with  the  corresponding  linear  statistical  model 

=  E  HiXU  +  eij 

e=i 

where  the  transformed  dependent  variables  are 


and  e jj  is  the  deviation  for  the  i^  sample. 
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The  strength  of  each  isotope  for  the  sample  is  estimated  by  the 
weighted  least-squares  method,  where  each  weight  is  inversely  proportional 
to  the  variance  of  the  datum.  Since  the  variances  are  unknown,  they  must 
be  estimated  from  past  and  current  data.  This  leads  to  a  nonlinear  problem 
which  is  solved  iteratively  by  starting  with  some  initial  estimates  for  the 
weights. 

Using  estimated  weights,  for  each  i  at  each  iteration,  minimize 
J  i 

2  "ijkhij  -  ^.Ik  (tjj,  T,j) 

i-1  L 

with  respect  to  /££ .  to  obtain  estimates  /Igik  by  solving  the  resulting  L 
simultaneous  linear  equations. 

The  estimation  of  weights  at  each  iteration  requires  the  partitioning  of 
the  estimated  variance  into  its  components.  As  a  consequence,  a  new  estimate 
of  the  sampling  error  os  is  obtained. 

An  estimate  of  the  total  variance  due  to  both  mechanical  and  Poisson 
sources  for  observation  of  the  ith  sample  is  the  usual 


'ijk  i?ij 


E 


namely,  the  error  variance  stemming  from  Poisson  phenomena  estimated  to  be 


plus  the  error  variance  stemming  from  mechanical  causes  estimated  to  be 


A  point  by  point  partitioning  of  each  deviate  into  amounts  ascribable  to 
mechanical  and  to  Poisson  sources  is  performed  to  obtain  a  point  estimate 
of  the  current  sampling  error  om. .  These  point  estimates  are  averaged  over 
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all  current  observations  to  get  am*k2,  which,  together  with  the  historic  estimate 
omo2,  are  averaged  to  provide  an  updated  historic  estimate  for  the  next  iteration. 
The  foregoing  procedures,  as  well  as  those  for  estimating  the  weights  and  source 
strengths,  are  summarized  in  the  following  simultaneous  equations: 

For  k  =  0, 1,2, ...,  last 


Km 


W 


ijk  +  l  '  N, 


Nil 


for  k  =  0 


■y—  +  T7  + 

1  ij  1 1| 


(1) 


K, 


ik  +  1 


Nii+M..k«tljfTii»  .  N3,C._] 

- r -  +  TT  +  °mk/i..k  (tjj,  Tjj ) 

1  ij  1 H  L  J 


for  k  >  0 


where  Kik  is  consistent  with 


J| 

E  wijk  + 1 

i-  1 


=  1 


(2) 


In  matrix  notation,  for  h  =  1,2, ...,  L, 


^Eik  +  1^  =  /l)  wijk+1Xhij  Xgij\  LE  wijk+1  X£ij  ^jj 


Lx  1 


'm'k  +  1 


2  _ 


E  <Ji  -  L> 

i= » 


j,  A 


-L  E 


ijk  +  1 


^ij  ^.ik+  1  ^ij» 


12 


L  ^.ik+  1  (tjj,  Tij)  J 


ft 


(3) 


E  <Ji  -  L> 

i  =  1 


where  ft  = 


[omk?.ik+l(tij,Tij)] 


°mk  M.ik  +  1 


2  M.ik  +  1  <^ij«  Tjj)  +  N*j 
Tii 


T*  U  • 

ij  wij 


12 


(4) 


*  ~  2  i  i  ~  2 

~  2  _  To  ®mo  +  T1  ®m‘k+ 1 

®mk  +  1  —  7  .  x 

To  T  T1 

where  amo  is  the  historic  estimate  of  the  proportional  mechanical  error  based 
on  f0  degrees  of  freedom,  and 


fi 


i 


E 


(Ji  -  L) 


is  the  additional  degrees  of  freedom  available  from  the  current  data. 

The  introduction  of  the  integers  Ujj  associated  with  background  reading 
requires  explanation.  A  single  background  reading  may  at  times  be  used  to 
correct  more  than  one  observation  of  the  same  sample  or  even  of  different 
samples.  The  product  of  the  number  of  samples  for  which  a  background 
reading  has  been  used  to  compute  tty's  times  the  number  of  observations  of 
the  ith  sample  for  which  this  particular  background  reading  has  been  used  is 
designated  Ujj.  It  is  introduced  into  the  iteration  merely  to  assure  that  the 
variance  due  to  a  single  background  reading  gets  counted  only  once  in  the 
averaging  process  over  all  observations  for  all  samples. 

Each  deviate  was  divided  by  /T  ik  (tjj(  Tjj)  to  assure  that  om.  will  be  on  a 
proportional  mechanical  error  basis,  as  defined.  The  last  multiplier  in 
Equation  3  partitions  each  squared  deviate  into  the  fraction  blamed  on 
mechanical  error.  The  sum  of  the  Jf  -  L  (degrees  of  freedom  associated 
with  ith  sample)  appearing  in  the  denominator  and  the  multipliers  Jj  -  L  in 
the  numerator  average  om2  over  all  the  total 

E  <Ji  -  L> 

i=  1 


degrees  of  freedom  available. 

Iteration  is  terminated  when  two  successive  am.  differ  by  less  than  an 
arbitrarily  small  amount,  6,  say  0.00005  (0.005%);  that  is, 

|  am'k  +  1  -  am'k  |  <  ^ 

At  the  end  of  the  least-square  fitting  process,  in  addition  to  the  final  am  , 
there  are  available  for  each  sample  the  final  L  +  1  estimates  contained  in 

-  E 


These  pg  (  and  I  final  inverse  matrices  with  associated  S(2  from  the 
least-square  fitting  process  are  used  in  establishing  confidence  intervals  and 
in  estimating  at. 

Estimation  of  p°  and  a, 

The  I  samples  have  been  used  to  obtain  the  final  estimates  p® ,  p  2< 
jTJ,  together  with  the  l-associated  final  L  by  L  inverse  matrices  and  Sj2  used 
for  the  least-squares  fits. 

In  using  these  it  will  be  assumed  erroneously  that  the  estimates  p0, , 
p^  , ....  p®  are  independent.  The  slight  dependency  that  may  exist  between 
two  particular  p°' s  stems  solely  from  the  use  of  some  of  the  same  background 
readings  for  the  computation  of  each  of  the  two  p“ 's.  This  dependency  is 
negligible  in  general  because  almost  invariably  N*  <  Njj.  To  consider  the 
dependency  would  unduly  complicate  the  computation. 

From  the  last  iteration  of  the  preceding  section 


where  Aj  is  the  sum  of  all  the  L°  elements  of  the  last  inverse  matrix  for  sample  i. 

To  obtain  p°  and  p°  oJ(  weights  v, ,  v2, v,  will  be  iterated  with  the 
l-fixedS^..2  and  'he  l-fixed  p°  as  inputs  to  obtain  a  stable  estimate  of  a.-.  The 

p.i 

iteration  algorithm  to  accomplish  this  consists  of  the  successive  use  of  the 
following  simultaneous  equations.  For  k  =  1 , 2, 3, ....  last 


v  =  -  for  k  =  1 

.-wo  \2  ,  o  2 

aso)  +  S_0 

M.i 


-  for  k  >  1 

/'■**-' O  ~  1 2  i  r»  2 

^..k-1  ask-l)  + 

H.\ 


where  Kk  is  consistent  with 
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Vk 


2 


/  ~  I4 

(K.kO(k-,) 


(?.°kO«k-i)  +S.02 

M.i 


(7) 


(I  -  1)os.k2  +  d  aso2 
d  +  I  -  1 


(8) 


where  aso  is  the  historic  estimate  of  proportional  samplir  g  error  based  on  d 
degrees  of  freedom.  The  historic  estimate  '.terns  from  previous  analyses  of 
all  pertinent  counting  data  from  a  laboratory  before  the 

E  Ji 

i*  i 


observations  currently  being  analyzed. 

Estimation  of  /ig 

Inputs  to  the  estimation  of  /ig 

1.  v, ,  v2, v,  from  the  last  iteration  of  the  preceding  section. 

2.  The/Tgj  (L  I  of  these)  from  the  final  iteration  of  the  section 
entitled  "Estimation  of  |ig(  and  om," 

Each  of  the  /To  will  be  computed  as  a  simple  partitioning  of  the  l i°  so 
that 

L 

~0  '^■'O 

-  E 

£-i 

I 

as  follows:  =  E  vi  Pg i 

i*i 


This  partitioning  is  based  not  on  statistical  grounds,  but  on  the  physical  fact 
that  all  samples  of  equilibrium  solutions  of  mixtures  contain  the  same  ratio  of 
ingredients,  to  a  precision  on  the  order  of  the  reciprocal  of  the  square  root  of 
the  number  of  molecules  involved  in  the  ingredient  sample. 
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Confidence  Intervals  for  n°  and  /ig 


Let  t*  be  a  value  computed  such  that 


Ut)dt 


0.95 


q  where  f(f)  is  the  Student's  t  density  function  for  t  with  I  -  1  degrees  of  freedom. 
Then  the  0.95  confidence  interval  for/i0  extends  over  the  interval 


0*0 


± 


and  for  jig  over  the  interval 


"C. 


± 


Partitioning  o~02 

The  estimated  variance  in  estimating  jj°  is 


~  2 
O^o 

M.. 


which  has  three  components  that  are  assumed  independent.  These  components 
are  errors  attributed  to  sampling,  mechanical  difficulties,  and  Poisson  phenom¬ 
ena,  and  may  be  expressed  as 


(v). + (v)„  +  (v), 


The  computation  of  these  three  components  will  now  be  discussed.  A 
crude  estimate  of 
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(>)  *  M 


E  =  £  v^AjS,2 


where  Af  is  the  sum  of  all  the  L2  elements  in  the  last 
and 

c2  1  i  r 


inverse  matrix  for  sample  i 


jpri  E  w‘i J^‘i  - 


To  partition  Sj2  into  its  Poisson  and  mechanical  components  and 
simultaneously  assure  that  the  sum  of  the  parts  equals  the  whole,  it  is 
necessary  to  partition  the  individual  squared  deviates.  The  point’ estimate 

for  the  jth  observation  on  the  i»h  samp|e  for  the  variation  due  to  the  Poisson 
phenomenon  is 

_  M.i  I«ij.  T,j)  +  Nil  n; 

ii  ^  ~ — 


uii  Tij 


and  for  the  variation  due  to  mechanical  causes 


fci  =  [®m?.i(tjj.TiJ)J2 

(°>°)  ’  ^  ^ 


which  guarantees  that 


Hence,  define 


where  a  and  p  are  crude  estimates  of 


,  respectively. 


~  2 

There  is  no  guarantee  that  a  +  p  may  not  exceed  o^0  ,  because  through  the 


vagaries  of  chance  the  I  values  of  y A  may  be  unusually  close  together. 


Finally 


a 

a  +  P"  +  7.. 


-  2 
0 — o 

M.. 


(9) 


P„ 

a.  +  P  +7., 


-  2 
O' — o 

V.. 


7..  _  2 

a  +  (3  +  7  M.. 


(10) 

(11) 


Summary 

The  Problem  Solution  section  may  be  summarized  briefly  in  the  following 
manner: 

1 .  T ransform  the  raw  data  as  specified  in  the  section  entitled, 
"Estimation  of /Ltg  and  am.” 
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2.  For  each  i  compute  a  least-square  fit  of  the  form 


Tjj)  =  Mli  +  ^2i  ^2ij  +  •••  +  ?Li  Xjj 

3.  From  these  I  results  estimate  om?  by  using  Equation  3. 

4.  Use  Equation  4  to  compute  om2  and  then  an  array  of  new 
weights  wjj2. 

5.  Continue  this  2-3-4  iterative  sequence  until 

|  k  |  <  5  x  fO'5 

6.  From  the  last  iteration  obtain  om. ,  am  ,  f0  +  fi ,  all  jig . ,  all  , 

and  all  S^,0  2 . 

M.i 


7.  Using  the  iterative  scheme  given  by  Equations  5  through  8, 
obtain  ?°  ,  as- ,  as ,  and  d  +  I  -  1  on  convergence;  that  is,  such  that 

|  <Vk  -  osk-i|  <  5x  10-5 

8.  Using  the  weights,  v/s,  of  the  last  iteration  in  step  7,  compute 
all  with  0.95  confidence  limits  on  JH°  and  each  /T°g  . 

9.  Use  Equations  9,  10,  1 1  to  partition  the  estimated  variance  of 
/j°  into  the  components  attributed  to  Poisson,  mechanical,  and  sampling 
causes. 


DISCUSSION  OF  MATHEMATICAL  MODEL 

The  procedure  outlined  in  this  report  estimates,  simultaneously,  total 
random-and-systematic  error,  its  Poisson  component,  its  sampling  component, 
and  its  residual  component,  called  mechanical  error  (a  presumed  composite 
from  all  causes  other  than  sampling  variation  and  Poisson  phenomenon).  This 
analytical  procedure  is  expressed  in  terms  of  the  error  components  at  time  t^, 
that  is,  age  zero.  The  error  under  discussion  is  error  in  estimating  corrected 
system  counts  per  minute  for  an  exact  intended  aliquot. 


In  principle,  this  unified  simultaneous  approach  shou\.  be  superior  to 
separate  estimations  of  the  three  foregoing  components  of  error.  In  passing, 
it  should  be  noted  clearly  that  this  report  is  concerned  with  an  analysis,  not 
a  calibration;  that  is,  laboratory  biases  will  never  be  brought  to  light  by  this 
procedure  alone.  The  analysis  may,  however,  be  used  for  the  calibration  of 
a  laboratory  when  laboratory  data  from  a  known  standard  solution  mixture 
are  available. 

Because  the  procedure  developed  in  this  report  is  an  initial  effort,  it  is 
not  the  last  word.  An  effort  was  made  to  obtain  either  unbiased  estimates, 
maximum  likelihood  estimates,  or  both,  so  that  the  results  could  be  claimed 
to  be  "best"  in  some  sense.  Unfortunately,  such  an  effort  did  not  seem 
feasible  at  all  points  of  the  development.  The  slight  statistical  imperfections 
that  have  crept  into  the  work  will  now  be  discussed. 

The  mathematical  model  assumes  that  the  standard  deviation  in  counts 
per  minute  for  the  random  error  from  all  causes  other  than  sampling  and 
Poisson  phenomenon  is  proportional  strictly  to  expected  corrected  counts 
per  minute  regardless  of  age,  isotope  solution  component  strengths,  or  length 
of  observing  period.  Perhaps  this  is  not  true  in  practice.  All  other  unaccounted- 
for  errors  are  encompassed  by  am  ,  so  that  (n°  om)2  plays  the  role  of  a  residual. 
Nevertheless,  the  unitless  proportionality  constant  om  is  very  meaningful  for 
the  physicist.  The  statistician  will  note  that  om  is  the  coefficient  of  variation 
for  "other  causes,"  assumed,  perhaps  erroneously,  to  be  constant. 

On  the  other  hand,  the  model  n°  os  for  standard  deviation  of  sampling 
random  error  should  be  almost  true  in  fact. 

In  computing  the  variance  of  a  datum,  expressed  as  counts-per-minute 
squared,  heavy  use  has  been  made  of  the  fact  that  for  a  Poisson  density  of 
counts,  the  mean  and  variance  are  equal.  Statisticians  generally  accept  this 
"fact,"  but  the  physicist  will  note  immediately  that  the  units  are  wrong; 
counts  cannot  equal  the  square  of  the  number  of  counts.  To  make  physical 
sense,  the  statisticians'  equation  needs  a  unit,  u*,  with  dimension  reciprocal- 
counts  so  that 


which  translated  says,  "for  the  population  of  a  fixed  Poisson  process,  the 
positive  square  root  of  the  numerical  value  of  the  mean  expressed  as  a  pure 
number  is  the  same  as  the  standard  deviation  of  the  population  expressed  as 
a  pure  number."  The  unit  u*  has  been  deliberately  left  out  of  the  Problem 
Solution  section.  To  check  on  unit  consistency  of  the  equations,  u*  will  need 
to  be  introduced. 
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The  Student's  t  density  has  been  freely  used  in  computing  confidence 
intervals,  even  though  it  is  realized  that  the  uncertainty  in  estimating  the 
various /i°'s  is  not  strictly  represented  by  a  normal  distribution,  since  some 
of  this  uncertainty  stems  from  nonnormal  Poisson  processes.  It  is  believed 
that  use  of  the  Student's  t  does  not  result  in  very  serious  errors.  This  belief 
is  based  upon  the  rapid  convergence  usually  encountered  in  applying  the 
central  limit  theorem  and  upon  the  fact  that  the  /j's,  J;'s,  and  particularly 

i 

Y  Jj  are  sufficiently  large  to  nearly  normalize  the  uncertainty. 

i  =  1 

Criticism  could  be  directed  at  the  following  type  of  sequential  steps 
which  weie  employed  freely: 

1.  oa+b2  =  aa2  +  ob2 

2.  oa+b2  =  oa2  +  ob2 

In  the  first  place,  1  is  not  true  unless  the  random  variables  a  and  b  are 
independent  of  each  other.  Where  1  was  used  the  two  random  variables 
involved  were  only  very  nearly  independent.  The  step  1  to  2  is  not  valid  in 
general  unless  the  quantities  in  2  are  unbiased. 

Additional  criticism  might  be  directed  at  partitioning  for  om  by  ratios 
rather  than  differences;  partitioning  os  is,  of  course,  open  to  the  same  criticism. 
This  method  of  partitioning  was  employed  to  avoid  negative  om2  and  negative 
os2,  since  such  results  are  both  impossible  and  meaningless  in  practice  and 
cause  insurmountable  difficulties  in  using  weights  based  on  estimated  variances. 
The  extent  of  possible  bias  introauced  by  this  method  cf  partitioning  is  un¬ 
known.  An  adequate  treatment  of  estimated  weights  or  of  partitioning 
variances  into  their  various  components  by  the  ratio  method  is  not  given  in 
the  literature,  and  investigation  of  these  subjects  would  constitute  a  larger 
study. 

In  spite  of  the  foregoing  reservations  it  is  believed  that  this  combined 
unified  approach  will  prove  to  be  a  significant  improvement  over  piecemeal 
analyses,  and  a  vast  improvement  over  assuming  that  all  variation  is  Poisson. 

Use  of  the  combined  approach  will  certainly  enable  a  laboratory  to  keep 
records  of  its  as  and  om  estimates  with  the  hope  that  these  estimates  will 
drift  lower  with  the  age  of  the  laboratory  as  techniques  improve.  Certainly 
these  estimates  will  point  to  the  area  where  the  most  improvement  is  to  be 
made.  Should  a and  am'  regress  significantly  with  time,  their  regression 
estimates  would  be  the  ones  to  use  in  the  weighting  process. 
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It  is  noted  that  under  many  sets  of  conditions  even  the  best  estimate  has 
an  appreciable  chance  of  missing  considerably.  For  the  problem  of  this  report, 
the  following  will  tend  to  improve  the  estimate  of  expected  decay  rate  at  some 
historic  time  of  interest: 

1  Decrease  of  sampling  error. 

2.  Sampling  replication. 

3.  Decrease  of  other  errors. 

4.  Replicating  and  lengthening  periods  of  observing  samples. 

Most  notable  is  item  2.  Presuming  reasonable  efforts  in  areas  1  and  3  have 
been  made,  but  an  appreciable  sampling  error  persists,  the  greatest  gain  per 
unit  of  effort  likely  can  be  accomplished  by  increasing  the  number  of  samples 
of  the  isotope  solution.  Simply  replicating  the  number  of  observations  of  a 
given  sample  will  not  average  out  sampling  fluctuation.  When  records  of  om 
and  as  are  in  hand  it  will  be  possible  to  use  them  to  design  a  schedule  for 
minimum-variance  estimate  within  a  fixed  total  cost. 


SIMULATION 

A  FORTRAN  program  was  written  to  implement  the  procedure  given  in 
the  Problem  Solution  section.  It  is  presented  and  discussed  in  the  Appendix. 
How  well  does  this  procedure  work?  In  order  to  obtain  a  preliminary  answer 
to  this,  it  is  desirable  to  operate  the  program  using  input  data  from  a  decay 
phenomenon  with  known  parameters.  To  accomplish  this  it  was  necessary  to 
develop  a  data  generator. 

Data  Generator 

The  known  constants  and  parameters  required  for  generating  data  are 
Xp  ,  L, /ig  ,  os,  crm,  I,  Jj,  tjj,  T|j,  T|*,  and  B (expected  background  counts  per 
minute). 

The  objective  is  to  simulate  I  samples  and  create  Jj  observations  (i  =  1 , 

2, ....  I)  at  different  ages  of  the  samples,  assuming  a  constant  parameter  for 
the  Poisson-distributed  background,  and  also  to  create  an  independent  set  of 
background  observations.  This  is  accomplished  in  the  following  steps. 

Step  1 :  I  random  numbers  are  drawn  from  a  normal  population  with 
mean  zero  and  standard  deviation  os.  Name  these  r, ,  r2, r, .  For  each  i 
compute: 
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IK7!  ■HMHKI 


t — "'"■■’■T.jr.jrrT 


Mu  -  <1  +  r.)  #«• 
A*2i  =  ^  +  ri)  M2. 

Mu  =  <1  +  M/Il. 


Step  2:  For  each  of  the  I  samples  draw  Jj  random  numbers  from  a 
normal  population  with  mean  zero  and  standard  deviation  om  .  Name  these 

1 

Qii »  Qi2  <  — «  Qijj .  Compute  ^  Jf  quantities 

i-  1 


£ij  “  (1  +  Qjj ) 


*xi  *'>  + 


MLi 


1-e'XLlii  \e-K^ 

VTij 


Step  3:  With  each  of  the  T jj  (£jj  +  Bjj)  as  the  Poisson  parameter,  select  at 
random  one  count  and  name  it  Cjj.  Select  at  random  one  count  from  each 
Poisson  density  with  parameter  T*  and  name  it  C(*. 

Fifteen  sets  of  data  were  generated.  For  each  set  the  input  constants  and 
parameters  were 


X,  =  0.30/min 
X2  =  0.15/min 
X3  =  0.05/min 
fi°  =  8,000  cpm 
H2,  =  1,500  cpm 
=  500  cpm 
as  =  0.05 

°m  -  0.08 

I  =  5 

Ji  =  10 

t jj  =  (0,2, 4, 6, 8, 15, 25,  35,45, 55)  minutes  every  i 
Tij  =  (1, 1, 1, 1, 1, 5, 5,  5,  5, 5)  minutes  every  i 
Tj|  =  10  minutes  every  i  and  every  j 


''l 

} 
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Bjj  =  10  cpm  every  i  and  every  j 
Ujj  =  1  every  i  and  every  j 

The  output  Cjj  and  Cj*  were  analyzed  with  the  computer  program  SAND 
(Appendix),  as  explained  in  the  next  subsection. 

Results 

The  15  sets  of  CSj  and  Cj*  together  with  the  information 
(\u\2,\3)  =  (0.30,0.15,  0.05)/min 
tjj  =  (0,  2,  4,  6,  8, 15,  25,  35,  45,  55)  minutes 
Tjj  =  (1,  1, 1, 1, 1,  5,  5,  5,  5,  5)  minutes 
Tj*  =  10  minutes 
t  =  0  minute 

were  inputted  to  the  computer  program  in  succession,  thereby  accumulating 
historic  am  and  os  with  associated  degrees  of  freedom  after  each  set  for  use  in 
the  next  set.  Set  1  was  rerun  after  set  15  and  then  set  2  was  rerun.  These 
runs  are  designated  V  and  2'  in  Tables  1  through  6,  and  the  slight  differences 
from  1  and  2  noted  stem  from  the  change  in  historic  om  and  as.  The  historic 
estimates  of  percentage  errors  are  graphed  in  Figure  1. 


Figure  1.  Estimated  sampling  and  mechanical  errors  versus  computer  run  number. 
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The  total  uncertainty  in  y°  expressed  as  a  variance  appears  in  the  last 
column  of  Table  1.  The  other  columns  contain  the  analysis  of  the  total 
uncertainty  into  sampling,  mechanical,  and  Poisson  sources. 

The  1 7  estimates  of  y° .  ,  and  with  0.95  confidence  intervals 

are  contained  in  Table  2.  It  will  be  recalled  that  y°  ~  10,000,  y°  =  8,000, 
y%  =  1,500,  and  y%  =  500  cpm. 

The  running  estimates  (current  and  historical)  of  om  and  os  are  contained 
in  Table  3.  The  degrees  of  freedom  for  1'  and  2'  are  admittedly  wrong  since 
the  old  data  of  1  and  2  were  used,  giving  really  no  new  degrees  of  freedom. 

In  detail  Table  4  shows  the  estimates  of  y°  ,  y°  ,  y^  ,  and  y\  from  each 
of  the  75  samples.  The  reader  will  note  the  entries  are  wilder  than  those  of 
Table  2,  which  showed  the  more  stable  weighted  averages  of  five  samples. 

Uncertainty  in  estimating  y°  from  a  single  sample  expressed  as  total 
variance  is  analyzed  into  mechanical  and  Poisson  sources  in  Table  5  for  each 
of  the  75  samples. 

For  each  of  the  75  samples  Table  6  displays  S|2  and  the  product  S;2  times 
the  elements  of  the  symmetric  inverse  matrix  under  the  headings  of  variances 
and  covariances. 


Table  1.  Analysis  of  Variance  of  Mean  Count 
Rate  at  Time  t  =  0  (All  Isotopes) 


Run 

Component 

Total 

Sampling 

Mechanical 

Poisson 

1 

234492. 

26212. 

14726. 

275432. 

2 

151186. 

34926. 

15107. 

201221. 

3 

6701 5. 

22903. 

8934. 

98853. 

4 

15476. 

7010. 

1553. 

24041. 

5 

6284. 

2416. 

646. 

9346. 

6 

43984. 

20101. 

6206. 

70292. 

7 

67485. 

37719. 

8886. 

114091. 

8 

148952. 

70537. 

16551. 

236041. 

9 

38675. 

14239. 

4133. 

57048. 

10 

566514. 

291941. 

96929. 

955385. 

11 

36820. 

19376. 

4379. 

60576. 

12 

136088. 

43760. 

11596. 

191446. 

13 

48631. 

27593. 

8649. 

84874. 

14 

23217. 

8967. 

1744. 

33929. 

15 

29841. 

15396. 

5465. 

50703. 

V 

180106. 

68296. 

11192. 

259595. 

2’ 

119166. 

65858. 

19668. 

204693. 

Table  2.  Mean  Decay  Rate  at  Time  t  =  0 
and  0.95  Confidence  Limits 


All  Isotopes  (cpm) 


Hun 

Lower  Limit 

Mean  (jj°) 

Upper  Limit 

9277.1 

10734.2 

12191.3 

8476.0 

9721.4 

10966.9 

8910.6 

9783.6 

10656.5 

9981.6 

10412.1 

10842.6 

9588.2 

9856.6 

10125.0 

9884.8 

10620.9 

11357.0 

8976.2 

9914.0 

10851.8 

8 

8759.3 

10108.2 

11457.0 

9 

9078.5 

9741.6 

10404.8 

10 

6461.8 

9175.6 

11889.4 

11 

9106.9 

9790.2 

10473.6 

12 

9254.6 

10469.4 

11684.2 

13 

9382.7 

10191.6 

11000.5 

14 

9662.2 

10173.6 

10685.0 

15 

8769.9 

9395.1 

10020.3 

V 

9461.6 

10876.2 

12290.8 

2 

8482.7 

9738.9 

10995.0 

First  Isotope  (cpm) 

Run 

Lower  Limit 

Mean  (p°) 

Upper  Limit 

7029.8 

8442.4 

9855.1 

6192.8 

7516.3 

8839.9 

6338.5 

752S  1 

8719.6 

7908.2 

8749.8 

9591.4 

7543.6 

7918.1 

8292.7 

7798.9 

8872.1 

9945.3 

7330.7 

7805.8 

8280.9 

8 

6886.9 

8200.7 

9514.5 

9 

6865.7 

7474.0 

8082.4 

10 

2678,6 

7037.2 

11395.8 

11 

6557.0 

7564.6 

8572.2 

12 

6460.0 

8587.5 

10715.0 

13 

7502.9 

8204.4 

8905.8 

14 

7036.4 

8069.6 

9102.7 

15 

6519.0 

7394.9 

8270.8 

V 

7177.7 

8577.6 

9977.4 

2 

6204.8 

7538.0 

8871.1 

Continued 
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Table  2.  Continued 


Second  Isotope  (cpm) 

Lower  Limit 

Mean  (ji^. ) 

Upper  Limit 

1663.5 

1914.2 

1180.2 

2279.1 

1203.0 

1754.9 

2306.8 

397.3 

1158.1 

1918.8 

1248.1 

wmm 

1616.9 

621.6 

1822.6 

906.2 

mm 

2333.5 

984.4 

1424.4 

1864.5 

1677.6 

1794.2 

1910.9 

-110.2 

1632.5 

3375.3 

1225.2 

1718.5 

2211.7 

238.7 

1 354.0 

2469.3 

1241.7 

1495.1 

1748.6 

713.4 

1613.9 

2514.4 

1045.0 

1524.1 

2003.3 

1671.7 

1791.7 

1911.7 

1183.0 

1725.0 

2267.0 

Third  Isotofx;  (cpm) 

Lower  Limit 

Mean  (JLI3  ) 

Upper  Limit 

476.4 

502.8 

529.1 

437.6 

475.3 

513.1 

456.6 

499.5 

542.4 

463.3 

504.1 

544.8 

466.6 

505.9 

545.1 

481.0 

526.6 

572.3 

459.1 

488.3 

517.4 

425.1 

482.9 

540.8 

391.5 

473.3 

555.0 

437.4 

505.8 

574.3 

472.1 

507.1 

542.0 

457.5 

527.8 

598.2 

466.8 

492.0 

517.2 

477.9 

490.0 

502.2 

408.8 

475.9 

543.1 

471.4 

506.8 

542.2 

439.3 

475.8 

512.3 

Table  3.  Estimates  of  Current  and  Historic  Mechanical  and  Sampling  Errors 


Run 

Mechanical  Error 

Sampling  Error 

Current  (%) 

Historic  (%) 

D.F. -^2 

Current  (%) 

Historic  (%) 

D.F.l^ 

0 

_ 

8.00 

0 

_ 

5.00 

0 

2.84 

2.84 

35 

9.90 

9.90 

4 

7.21 

5.48 

70 

8.95 

9.44 

8 

9.81 

7.22 

105 

5.47 

8.33 

12 

5.82 

6.90 

140 

2.64 

7.33 

16 

6.09 

6.74 

175 

1.78 

6.60 

20 

7.78 

6.93 

210 

4.47 

6.30 

24 

6.28 

6.84 

245 

5.88 

6.24 

28 

8 

5.58 

6.69 

280 

8.44 

6.56 

32 

9 

8.03 

6.86 

315 

3.73 

6.31 

36 

10 

12.03 

7.53 

350 

12.16 

7.11 

40 

11 

8.13 

7.59 

385 

4.42 

6.91 

44 

12 

7.23 

7.56 

420 

7.95 

7.00 

48 

13 

9.63 

7.74 

455 

4.93 

6.87 

52 

14 

4.84 

7.57 

490 

3.34 

6.68 

56 

15 

9.44 

7.71 

525 

4.11 

6.54 

60 

r 

4.88 

7.56 

560 

8.49 

6.68 

64 

CM 

8.40 

7.61 

595 

7.96 

6.76 

68 

\J  D.F.  =  Degrees  of  freedom  associated  with  updated  historic  estimate. 


Table  4.  Mean  Decay  Rates  at  Time  t  =  0 


Run 

Sample  (i) 

Isotope  (cpm) 

All  (p°) 

First  (£?  ) 

Second  (pt^  ) 

Third  ) 

1 

1 

9075.5 

6698.7 

1859.0 

517.7 

2 

12219.9 

9790.3 

1901.1 

528.4 

3 

11225.4 

8920.0 

1813.9 

491.4 

4 

10018.4 

7880.4 

1663.9 

473.9 

5 

10774.7 

8572.0 

1701.5 

501.1 

2 

1 

9152.7 

7013.6 

1637.8 

501.2 

2 

10766.9 

8815.5 

’465.9 

485.4 

3 

10830.6 

8453.6 

1889,0 

487.9 

4 

8727.2 

7021.8 

1222.0 

483.3 

5 

9110.2 

6327.0 

2358.3 

424.8 

3 

1 

9923.6 

8238.1 

1157.8 

527.6 

2 

9565.3 

6857.8 

2215.4 

492.0 

3 

11121.8 

9055.2 

1521.8 

544.7 

4 

9249.4 

7101.5 

1693.7 

454.1 

5 

9366.8 

6755.2 

2121.8 

489.7 

Continued 
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Table  4.  Continued 


Isotope  (cpm) 


Sample  (i) 


All  (£°) 

First  0T°) 

Second  (#/2. ) 

Third  (£3^ 

10597.1 

8923.6 

116'  2 

507.2 

10006  6 

7925.4 

1602.2 

478.9 

10483.6 

9737.8 

187.6 

558.0 

10137.5 

8493.7 

1154.0 

489.7 

10872.9 

8647.2 

1742.3 

483.3 

9694.7 

7756.2 

1440.8 

497.6 

10051.1 

8346.1 

1180.6 

524.3 

9980.2 

7846.5 

1583.2 

550.4 

9565.7 

7566.8 

1514.8 

483.9 

10019.6 

8130.0 

1417.2 

472.2 

10289.3 

8260.2 

1563.1 

465.8 

10728.6 

8634.0 

1518.6 

575.9 

9767.4 

8194.1 

1057.4 

515.8 

11293.6 

10248.6 

510.2 

534.7 

10852.0 

8724.1 

1596.6 

531.2 

9463.3 

7684.7 

1293.5 

484.9 

10570.3 

8039.8 

2026.3 

504.0 

9670.3 

7307.0 

1910.5 

452.8 

9013.5 

7779.7 

727.2 

506.6 

i  0775.2 

8249.1 

2026.4 

499.6 

8834.8 

7204.9 

1188.8 

440.9 

11592.1 

9376.7 

1688.5 

526.8 

9309.6 

7009.9 

1872.8 

426.8 

10418.7 

8956.2 

954.8 

507.6 

10558.5 

8676.3 

1365.4 

516.7 

9787.4 

7681.8 

1719.1 

386.4 

10883.7 

8397.0 

1990.4 

496.2 

9157.6 

6869.5 

1804.1 

483.8 

9876.2 

7476.4 

1854.9 

544.7 

9601.5 

7340.8 

1744.5 

516.2 

10068.6 

8732.8 

786.8 

548.9 

2386.6 

-4007.7 

6037.1 

357.2 

9539.8 

7265.5 

1781.0 

493.2 

8612.3 

6462.7 

1669.7 

479.7 

10428.2 

8743.2 

1147.4 

537.4 

9303.7 

7232.3 

1586.9 

484.5 

9911.3 

7996.5 

1411.6 

503.0 

10604.2 

8593.1 

1488.7 

522.3 

9275.1 

6340.7 

2454.2 

480.1 

9824.4 

7339.7 

1932.4 

552.2 
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Table  4.  Continued 


Run 

Sample  (i) 

Isotope  (cpm) 

AM  <M°) 

First  <^i1  ) 

Second  (*i2  ) 

Third  (pt j  ) 

12 

1 

12061.1 

11422.0 

23.0 

616.0 

2 

10242.7 

8069.6 

1663.2 

509.8 

3 

9720.2 

8114.4 

1080.6 

525.1 

4 

10308.4 

7497.9 

2349.3 

461.1 

5 

9761.9 

7316.9 

1936.7 

508.2 

13 

1 

9169.9 

7413.9 

1253.7 

502.2 

2 

10490.0 

8786.0 

1254.7 

449.1 

3 

10880.1 

8763.3 

1621.2 

495.5 

4 

10218.1 

8104.0 

1613.8 

500.2 

5 

10274.7 

8110.2 

1663.9 

500.5 

14 

1 

9786.2 

8217.4 

1074.2 

494.4 

2 

9683.3 

7503.2 

1700.3 

479.7 

3 

10645.6 

9381.8 

781.9 

481.8 

4 

10361.3 

7324.9 

2531.9 

504.3 

5 

10349.2 

7743.9 

2115.4 

489.9 

15 

1 

8783.8 

6940.2 

1377.3 

466.1 

2 

9021.0 

6618.1 

1957.9 

444.9 

3 

9862.1 

8028.2 

1282.4 

551.4 

4 

9866.4 

8366.9 

987.3 

512.0 

5 

9640.6 

7392.1 

1831.3 

417.1 

r 

i 

9326.7 

7014.7 

1785.1 

526.9 

2 

12396.6 

10078.2 

1775.5 

5^2.8 

3 

11184.0 

8774.0 

1927.7 

482.2 

4 

10081.1 

7804.9 

1801.4 

474.7 

5 

10853.6 

8677.5 

1672.6 

503.4 

2' 

1 

9159.6 

7017.2 

1641.4 

500.9 

2 

10791.1 

8870.5 

1431.1 

489.4 

3 

10835.2 

8457.6 

1890.0 

487.5 

4 

8727.7 

7007.6 

1238.0 

482.1 

5 

9129.1 

6370.9 

2330.7 

427.4 
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Tabic  5.  Analysis  of  Variance  of  Mean  Decay  Rate  at  Time  t  =  0 
(All  Isotopes  -  Fixed  Sample) 


Run 

Sample  (i) 

Weight  (Vj) 

Component  of  Variance 

Total  Variance 

Mechanical 

Poisson 

1 

1 

0.1784 

263289. 

103927, 

367217. 

2 

0.2106 

88950. 

49163. 

138114. 

3 

0.2142 

71187. 

45703. 

116891. 

4 

0.1829 

187432. 

142537. 

329970. 

5 

0.2136 

78093. 

42290. 

120384. 

2 

1 

0.1834 

234844. 

149464. 

384308. 

2 

0.1928 

233754. 

90652. 

324406. 

3 

0.2093 

181006. 

51120. 

232127. 

4 

C  .1990 

234582. 

53492. 

288075. 

5 

C.2153 

113077. 

89236. 

202314. 

3 

1 

0.2096 

212042. 

69042. 

281085. 

*7 

0.2214 

184160. 

46395. 

230555. 

3 

0.1599 

391622. 

183624. 

575247. 

4 

0.2091 

183676. 

99464. 

283140. 

5 

0.1997 

238789. 

89202. 

327992. 

4 

1 

0.2076 

231941. 

60566. 

292507. 

2 

0.2114 

254415. 

22040. 

276455. 

3 

0.2081 

243324. 

47046. 

290371. 

4 

0.1894 

320292. 

56082. 

376374. 

5 

0.1833 

286201. 

122176. 

408378. 

5 

1 

0.2254 

91493. 

47235. 

138729. 

2 

0.1823 

235998. 

35760. 

271759. 

3 

0.2068 

162977. 

26366. 

189344. 

4 

0.1924 

172516. 

62571. 

235088. 

5 

0.1929 

186934. 

46729. 

233664. 

6 

1 

0.1604 

331838. 

127761. 

459599. 

2 

0.1835 

297079. 

48237. 

345316. 

3 

0.2012 

180200. 

95300. 

275501. 

0.226"' 

151903. 

42326. 

194230. 

!  > 

0.r.279 

153969. 

36814. 

190784. 

7 

i 

0.1589 

387264. 

67945. 

455209. 

2 

0.1648 

317444. 

107897. 

425341. 

3 

0.2380 

159895. 

16860. 

176756. 

4 

0.2019 

200948. 

75684. 

276633. 

5 

_ 

0.2361 

150803. 

30359. 

_ 

181162. 

Continued 
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Table  5.  Continued 


Run 

Sample  (i) 

Weight  (Vj) 

Component  of  Variance 

Total  Variance 

Mechanical 

Poisson 

8 

1 

0.2135 

181379. 

45646. 

227026. 

2 

0.1840 

307707. 

26350. 

334058. 

3 

0.2005 

215719. 

54452. 

270172. 

4 

0.1569 

383884. 

83341. 

467225. 

5 

0.2448 

95522. 

46152. 

141674. 

9 

1 

0.2914 

34774. 

7777. 

42551. 

2 

0.1112 

549551. 

173906. 

723458. 

3 

0.2207 

146099. 

31088. 

177187. 

4 

0.1502 

294614. 

142978. 

437592. 

5 

0.2262 

137816. 

25790. 

163606. 

10 

1 

0.2698 

116786. 

54891. 

171678. 

2 

0.0682 

1 752976. 

184104. 

1937081. 

3 

0.1959 

354252. 

43133. 

397385. 

4 

0.2380 

162209. 

89331. 

251540. 

5 

0.2278 

189679. 

92461. 

282140. 

11 

1 

0.2535 

102776. 

53193. 

155969. 

2 

0.2177 

232965. 

23992. 

256957. 

3 

0.2088 

244970. 

42622. 

287593. 

4 

0.1527 

481337. 

80013. 

561350. 

5 

0.1670 

367288. 

106687. 

473975. 

12 

1 

0.2192 

149618. 

23907. 

173525. 

2 

0.2579 

47669. 

18974. 

66643. 

3 

0.1851 

253867. 

50431. 

304299. 

4 

0.1596 

321370. 

117524. 

438895. 

5 

0.1779 

265270. 

73032. 

338302. 

13 

1 

0.2122 

264342. 

64764. 

329106. 

2 

0.1501 

536774. 

131538. 

668313. 

3 

0.2172 

210208. 

100169. 

310378. 

4 

0.2188 

234266. 

70318. 

304584. 

5 

0.2014 

280819. 

92317. 

373136. 

14 

1 

0.2136 

136544. 

36041. 

172586. 

2 

0.1808 

265517. 

2241 1. 

287929. 

3 

0.2118 

131970. 

46039. 

178010. 

4 

0.1886 

223318. 

33582. 

256900. 

5 

0.2049 

167787. 

31691. 

199478. 
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Table  5.  Continued 


Run 

Sample  (i) 

Weight  ( V j ) 

Component  of  Variance 

Total  Variance 

Mechanical 

Poisson 

15 

1 

0.2183 

159831. 

63264. 

22309* 

2 

0.2219 

170052. 

43352. 

213405. 

3 

0.2123 

175151. 

64917. 

240068. 

4 

0.1417 

400372. 

147496. 

547869. 

5 

0.2057 

186634. 

73194. 

259828. 

r 

1 

0.1738 

336522. 

29850. 

366372. 

2 

0.2167 

159610. 

29845. 

189455. 

3 

0.2240 

131648. 

34580. 

166229. 

4 

0.1558 

393045. 

76615. 

469660. 

5 

0.2294 

134413. 

15233. 

149647. 

2' 

1 

0.1688 

331167. 

125479. 

456646. 

2 

0.1887 

287619. 

75461. 

363080. 

3 

0.2185 

206020. 

48241. 

254262. 

4 

0.2042 

265972. 

36648. 

302621. 

5 

0.2195 

160941. 

90226. 

251168. 
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Table  6.  Variances  and  Covariances  of  Fitted  Parameters  and  Weighted  Variance  of  Errors 
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Discussion  of  Results 


Fifteen  runs  ori  different  simulated  samples  of  data  with  parameters 
H  =  //,  +  n2.  +  M3,  equal  to 

10,000  =  8,000  +  1,500  +  500 

have  been  analyzed  by  means  of  the  FORTRAN  program.  For  each  run  0.95 
confidence  intervals  for  ,  n? .  and  /LI3  were  printed  as  outputs  (see 

Table  2).  If  the  procedure  were  perfect,  in  the  long  run  5%  of  the  intervals 
would  fail  to  bracket  (contain)  the  parameter  represented.  Five  percent  of 
fifteen  is  0.75.  Thus  a  priori  the  expected  number  of  misses  per  parameter 
is  0.75  and  for  all  four  parameters  three  misses.  Actually,  two  misses  were 
experienced,  ^  at  run  1  and  at run  9-  If  the  procedure  is  imperfect 
with  reference  to  confidence  intervals,  a  large  number  of  additional  runs 
would  be  needed  to  "prove"  it  statistically. 

The  data  generator  introduced  ...echanical  error  with  parameter 

am  -  0.08 

From  the  15  runs  with  their  available  525  degrees  of  freedom  for  estimating 
om  ,  the  final  estimate  was: 


om  =  0.0771 

Is  the  procedure  biased  downward?  The  difference  -0.0029  is  not  statistically 
significant  at  the  0.95  level.  Thus  more  runs  would  be  needed  to  prove  that 
the  procedure  is  imperfect  for  estimating  om  .  Even  if  the  procedure  is  imper¬ 
fect,  the  bias  appears  to  be  small.  The  authors  anticipated  a  small  upward  bias 
rather  than  the  downward  bias  actually  observed,  because  the  variance  of  the 
sum  of  the  mechanical  and  Poisson  components  should  be  slightly  more  than 
the  sum  of  the  two  components,  due  to  a  small  positive  covariance.  This 
covariance  stems  from  the  slight  tendency  for  large  Poisson  errors  to  be 
associated  with  large  mechanical  error  since  om  is  on  a  percentage  basis.  The 
procedure  ignores  this  tendency  by  partitioning  as  though  the  two  sources 
were  independent.  The  procedure  should  have  a  tendency  to  lead  to  a 
divisor  too  small  in  the  proportioning  for  om  and  thus  give  a  om  too  large. 

The  data  generator  introduced  sampling  error  with  parameter 

os  =  0.05 
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From  the  15  runs  with  their  available  60  degrees  of  freedom  for  estimating  as, 
the  final  estimate  as  =  0.0654  misses  the  parametric  value  by  a  disquieting 
amount.  While  the  physicist  may  not  be  concerned  with  a  bias  of  1-1/2%, 
this  is  statistically  significant  at  the  0.99  level,  from  which  it  is  concluded  that 
the  procedure  estimates  os  with  an  upward  bias.  Why?  Some  of  it  can  stem 
from  the  fact  that  the  procedure  treats  sampling,  Poisson,  and  mechanical 
error  sources  as  independent,  whereas  because  of  the  percentage  concept 
there  are  actually  three  small  positive  covariances  between  these  three  sources. 
In  the  proportioning  for  os  these  were  not  included  in  the  divisor  -  which 
should  result  in  a  tendency  to  biasos  upward. 

It  is  noted  that  os  is  estimated  at  the  last  step  of  the  computation.  It  is 
therefore  suspected  that  the  accumulated  effects  of  all  the  slight  imperfections 
of  the  procedure  finally  appear  in  os .  Two  of  the  approx  imations  suspected 
of  contributing  a  large  part  of  the  imperfections  are  ( 1 )  the  use  of  estimated 
variances  instead  of  true  variances  in  the  weighting  and  (2)  the  inefficient 
estimation  of  Poisson  parameters  by  the  partitioning  of  each  regression  error. 
(A  discussion  of  the  efficiency  of  estimating  Poisson  parameters  can  be  found 
in  another  report.*) 

Imprecise  as  the  procedure  may  be  it  is  believed  to  be  a  much  better 
approach  than  the  customary  assumption  that  all  error  is  Poisson  distributed. 

It  is  hoped  that  physicists  will  use  this  procedure  and  that  statisticians  will 
attempt  to  improve  it. 


FINDINGS 

A  procedure  implemented  as  the  SAND  computer  program  is  now 
available  for  use  by  nuclear  physicists.  Given  a  multiple  isotope  source  and 
a  series  of  observations  on  some  samples  of  the  source,  this  FORTRAN  com¬ 
puter  program  will  estimate  the  strength  of  each  isotope;  place  confidence 
limits  thereon;  break  the  uncertainty  into  Poisson  error,  sampling  error,  and 
residual  (called  mechanical)  error;  supply  an  estimate  of  the  standard  deviation 
of  sampling  error;  and  supply  an  estimate  of  the  standard  deviation  of 
mechanical  error. 


*  U.  S.  Naval  Civil  Engineering  Laboratory.  Technical  Report  R-519  Efficiency  of  two 
estimates  for  a  Poisson  distribution,  by  W.  L.  Wilcoxson  and  M.  L.  Eaton.  Port  Hueneme, 
Calif.,  Apr  1967. 
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CONCLUSIONS 


The  procedure  developed  in  this  report  constitutes  a  considerable 
improvement  over  the  customary  technique,  whereby  overoptimistic  con¬ 
fidence  intervals  have  been  based  on  only  Poisson  variation.  Even  though  the 
procedure  contains  statistical  imperfections  that  would  require  conside-  ibly 
more  research  to  remove,  it  does  a  reasonably  good  job  of  breaking  total 
variation  into  its  three  sources. 


RECOMMENDATION 

It  is  suggested  that  nuclear  physicists  use  the  procedure  described  in  this 
report  for  the  statistical  analysis  of  counting  data. 
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Appendix 


SAND,  A  FORTRAN  ll-D  COMPUTER  PROGRAM  FOR  THE 
STATISTICAL  ANALYSIS  OF  NUCLEAR  DATA 


The  procedure  described  within  the  body  of  this  report  is  implemented 
by  the  computer  program  SAND  on  an  IBM  1620  with  an  IBM  131 1  Disk  Drive. 
A  minimum  1620  (20K  memory)  with  no  special  features  but  a  disk  is 
required.  The  programming  language  is  the  disk  version  of  FORTRAN  II; 
however,  the  disk  requirement  may  be  replaced  by  magnetic  tape  by  changing 
the  FETCH  and  RECORD  statements  in  the  following  subroutines:  FANDW, 
INPUT2,  and  BCKGND. 

All  input  to  the  program  is  from  the  card  reader,  and  all  output  is  to  the 
card  punch  for  off-line  listing  on  an  IBM  407.  The  card  I/O  statements  may 
be  replaced  by  tape  I/O  through  a  simple  sifting  program. 

The  program  for  the  IBM  1620  is  limited  to  a  maximum  of  25  samples 
per  run. 

The  format  of  the  input  for  each  run,  the  card  deck  setup  (Figure  2),  a 
sample  input  with  corresponding  output  (Figure  3),  and  the  listing  of  the 
FORTRAN  I  l-D  source  program  (Figure  4}  follow. 


FORMAT  OF  INPUT  FOR  EACH  RUN 


Card  1.  Title  card;  FORMAT  (80H 

Card  2.  I,  L,  r,  tQ.  omo,  f0.  oso,  d,  amr ;  FORMAT  (215,  7E10.0) 


where  I  =  numbers  of  samples  (I  <  25) 

L  =  number  of  isotopes  (L  <  5) 
t  =  detector  deadtime  in  min/count 


tD  =  historic  time  at  t  =  0  (minutes) 

omo  =  historic  estimate  of  mechanical  error  (percent! 

fD  =  degrees  of  freedom  associated  with  omo.  ( I f  f 0  =  0,  then 
this  is  a  starting  estimate  used  in  solving  the  nonlinear 
equation  in  om.) 


I 

I 

I 


oso  -  historic  estimate  of  sampling  error  (percent) 

d  =  degrees  of  freedom  associated  with  oso.  (If  d  =  0,  then  this  is  a 
starting  estimate  used  in  solving  the  nonlinear  equation  in  os.) 

omr  =  partially  updated  omo  used  in  restarting  program  after  interruption 
(normally  blank) 

Card  3.  X,,X2 . XL;  FORMAT  (5E15.8) 

where  Xg  =  decay  constant  for  isotope  2  (2  =  1,2 . L  <  5)  in 

reciprocal  minutes 

Card  4.  J1(  J2 . J, ;  FORMAT  (1515) 

where  J(  =  number  of  observations  for  sample  i  (i  =  1,  2 . I  <  25) 

Observation  Cards.  Tjj,  Cjj,  Tt^,  C,* ,  ;  FORMAT  (5E15.8) 

where  the  subscript  j  varies  the  most  rapid;  j  =  1 ,  2,  3,  J  j  for  i  =  1 , 2, 


Figure  2.  Deck  setup  for  each  run. 
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Figure  3.  Sample  input  and  output. 
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Figure  3.  Continued. 
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Figure  3.  Continued. 
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Figure  3.  Continued. 
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Figure  3.  Continued. 
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Figure  4.  Program  listing. 
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Figure  4.  Continued. 
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Figure  4.  Continued. 
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Figure  4.  Continued. 
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Figure  4.  Continued. 


n 


c 


3 

a 

z 


<  — 


s 

•> 

-)  — 
—  cc 
5i  • 
>  iTi 

<  H 

•  Ul 
“>  <P 


<1 

D 


* 

H 

3 

a 

*— J 

h 

• 

U; 

J— • 

» 

2 

» 

« 

o 

LJ 

a 

3 

•> 

IT 

X 

Z 

— * 

CV. 

< 

Uj 

if\ 

w 

w 

> 

w 

X 

h- 

UJ 

•> 

X 

< 

♦ 

• 

c 

to 

•w. 

c 

2. 

u 

o 

K 

X 

•» 

z 

*) 

» 

3 

» 

u 

< 

QT 

x 

3 

UJ 

CO 

to 

< 

UJ 

< 

x 

u_ 

IT. 

K 

•» 

K 

Li_ 

c 

fVi 

10 

>- 

» 

» 

#k 

< 

w 

f- 

Q 

X 

to 

(- 

r 

3 

X 

a 

< 

£: 

a: 

s: 

x 

z 

o 

t- 

U 

< 

UJ 

3 

< 

X 

oa 

Uj 

»- 

* 

to 

tO 

a 

j- 

> 

lT 

i 

•> 

UJ 

< 

tO 

e 

u 

u 

CM 

» 

«— 

u, 

X 

#* 

U- 

tO 

U- 

c0 

_i 

UJ 

if. 

UJ 

z 

•> 

D 

L*J 

U' 

» 

f\> 

• 

Uj 

r-J 

>. 

» 

D 

to 

«■»-* 

*— « 

a 

* 

X 

3 

— * 

Z 

u. 

3 

K 

> 

X 

< 

to 

CO 

o 

Z 

♦ 

< 

o 

x 

» 

W 

u 

« 

u 

* 

#■ 

_i 

a 

X 

o 

to- 

l/. 

•* 

x 

c 

» 

z 

UJ 

► 

2 

> 

K- 

X 

c. 

#> 

to 

•> 

— 

UJ 

< 

_J 

<. 

to 

*— 4 

s: 

i r. 

•> 

»  UJ 

z 

u_ 

2; 

to 

Ll 

<► 

ir  .  a; 

> 

o 

_i 

O 

U- 

>- 

O 

ir, 

C\J  < 

#* 

•> 

#> 

#> 

Cj 

_J 

» 

c 

cc 

O 

u 

u. 

a 


m  “i-  <  a 

r  >  zr  z  < 

2  <1  LiJ  < 


•  > 

x  < 


X 

< 

51 

Vi 


—  IT 

\0  cm 


C  •  w 


a 

> 

< 


ir> 


a: 

O 


*  co  h- 


”>  — 


u 

UJ 


3 

a 

z 


«r  X 

<  w  > 
<  < 


cc 

O 


vO  IT 
«-  CM  < 


—  < 
a  x 

>  3  tO 

<  ►  o 

►  Da 
3  X 
<->  UJ 
X  tO  * 

>  >  51 
*  < 


X 

a 

< 

tO 

> 


CM 


< 

z 

< 


X1 

•* 

o 


z  -* 


X  — 
x  a  u  x 

<  Z  UJ  u 

Z  <  5  • 

>  ->  to  u.  — 

<  •>  >  O 

Q.  *  *  w 

Z  5  I  X 

<  X  X  — 


< 

u 


3 

< 


to 


< 
I- 
i O 


X 

< 

z 


e 

CM  > 

o 

10 

X 

w 

•* 

3  a 

lO  U 

u 

UJ 

r— 1 

< 

X 

UJ 

to 

z 

> 

5l 

> 

LJ 

UJ 

u 

I 

X 

X 

» 

X 

z 

X 

< 

z 

» 

UJ 

< 

r 

X 

X 

z 

in 

o 

(\J 

o 

a: 

o 

X 

to 

•  cc 

> 

> 

UJ 

O  4* 

o 

co 

if> 

W 

»- 

U  H 

J— 1 

IT 

» 

X 

•h 

3 

c 

CM 

— - 

r 

f-  Q 

z 

O 

c 

3 

3 

if) 

3 

u 

•* 

I 

Tj^ 

< 

(\J 

»- 

t— 

CO 

I 

< 

a: 

a: 

»- 

t- 

O 

*— t 

• 

z 

<  o 

X 

UJ 

X  u 

o 

> 

X  < 

X 

< 

U  X 

O 

3 

3 

3 

a 

Q 

u 

Q 

z 

UJ 

;r 

< 

UJ 

X 

*— 1 

DUX 

U  X 

Q 

z 

a 

u 

t- 

o 

CL 

a 

no 

< 

z 

2: 

X 

> 

V 

X 

3 

< 

z  cc 

z 

cc 

< 

3 

O 

UJ 

UJ 

z 

Z 

z 

3 

• 

UJ 

•— 1 

< 

o 

< 

< 

> 

o 

o 

UJ 

3  O 

3 

o 

UJ 

a 

u. 

O' 

or 

UJ 

M 

M 

to 

3 

cc 

c 

u 

u 

UJ 

a 

a  u. 

a 

u. 

a 

»— 

LC 

r-H 

H 

CM 

CO 

*H 

c 

o 

c 

o 

CC 

a 

UJ 

cc 

to 

o 

o 

IP 

o 

r~ 

o 

3 

-1 

o 

•— 

rH 

CM 

o 

Uj 

3 

Z 


Z 

o 

u 


ICUiU 
N  NQ  N 

M  K  *  M 


U  U 


54 


r-(\jrO.JiniONCDO-  C  fCMrO,*iniOr^OCO-CftCMm>Jin.OfvoO 
O'  O  O  O  O'  O  O  O  O  0000000000<-><»-'»-<*^*-*«—  hhh 

ftp—fti-ifififififi(Mr\fMtMr\j(MtMt\tCM<MCMrMCMfMfMCMCMCMCM 

cc  cccooocccccococcoccooccccc 


ro 


■g 


3 

c 


c 

o 

o 


'J- 

Q> 


3 

O) 


LO 

f 

U. 

L— 

a 

z 

f 

C 

X 

— ■ 

t- 

o 

CM 

o 

o 

• 

V 

vO 

•— 

z 

ft 

ft 

(J. 

Li. 

X 

*> 

f 

o 

X 

f 

I 

c- 

cc 

CL 

c 

r- 

ft 

U. 

ft 

w 

» 

(Li 

» 

V 

X 

X 

<M 

z 

lO 

X 

• 

ft 

CM 

z 

t4 

X 

«— 

ft 

Vs 

V 

Ll 

ft 

CO 

X 

« 

ft 

c. 

UJ 

I 

* 

CM 

X 

»o 

Ll 

V 

ft 

♦ 

_l 

«H 

1- 

CM 

a 

1“ 

1— 1 

c 

• 

X 

3 

*— 

< 

'O 

< 

a 

« 

IL 

ft 

lO 

z 

l/i 

Q 

a 

>— 

U- 

• 

CL 

-1 

CL 

tO 

■*-v 

U. 

a 

c 

O  - 

X 

1-  f 

X 

a 

3 

X 

II 

1—  X 

< 

Z  X 

< 

< 

< 

t- 

u  < 

X 

<  < 

10 

t- 

CO 

IL  X 

h-  X 

*-* 

z 

L— 

►_  _J 

#t 

co  •- 

♦ 

o 

c 

U. 

< 

lL  •> 

rH 

z  » 

ft 

1- 

O 

O  r-1 

II 

O 

II 

1- 

•t 

UJ 

II 

-1 

U  II 

•— i 

< 

— * 

X 

CL 

X 

_J 

• 

ft 

> 

c 

< 

Uj 

» — 1 

X  • 

>-  • 

a 

• 

J* 

a: 

1— 

O'  f 

-J 

<  - 

UJ 

c 

_i 

X 

ft _ i 

W 

L  J  *- 

w 

cO 

»H 

• 

D 

u 

V  - 

< 

UJ  f 

X 

cc 

L- 

X 

Z 

ft 

z  < 

CO  o 

o  X 

< 

o 

h- 

< 

CL 

—  D 

•  X 

< 

«"-v 

X 

•* 

X 

O  X  X 

in 

< 

X 

IT 

•J 

in 

1 — I 

X 

K 

< 

H 

w 

X  -> 

*— 1 

W 

i 

M 

•t 

O' 

CO 

X  ~ 

UJ 

* 

|V  — 

in 

O' 

LO 

r\j 

c 

»— H 

— 

<f  * 

i n  c 

ft  • 

«-t  o 

<\J 

W 

•J 

w» 

X 

•  c 

— 

»o 

w  C 

*— 

CO 

w 

H 

h- 

cm  m 

i- 

►—  r- 

1- 

1-  z 

< 

X 

< 

• 

< 

X 

< 

< 

X 

<  CL 

5. 

u 

X 

X 

rH  Q 

X 

u 

X  c 

X 

u 

X  X 

a 

z 

CL 

(M 

rH  <£ 

CL 

z 

CL  < 

CL  Z 

a  k 

o 

O  D 

O 

NU.LJ03 

O  UJ  o 

D 

O  UJ 

z 

U 

a 

Ll 

CL 

a 

0. 

U.  CL 

U. 

a 

a  a 

IL’ 

H 

<M 

C 

O 

o 

o 

O 

o 

a 

n 

■4- 

in 

vC 

|V 

CO 

o 

m 

cm 


CL 

X 

< 

lO 


X 

< 

X 

“> 


a 

< 


c 

D 

21 

IL.' 


D 

a 

z 


m 

K  • 

OJ 

CO  1- 

f 

1-  Q 

3 

•>  r> 

X 

a  x 

UJ 

<  a 

cO 

h-  * 

> 

LO  X 

• 

<sJ  U 

f 

•  a 

in 

a  x 

CM 

x  •  a 

— ■ 

ft  Q  «t 

i-  >  a 

X 

•  <  Q 

a 

D  *  » 

cO 

*  X  < 

« 

Z  >  1- 

a  <  _j 

m 

•  •  UJ 

UJ 

• 

in 

a  >  o 

z 

m 

CM 

4  »  » 

ft 

CM 

W 

i-  <  a 

h- 

f 

a 

to  o  X 

o 

D 

> 

Z  X  < 

C 

X 

< 

a  <  to 

a 

a 

* 

*  •  > 

• 

f 

->  x  < 

i- 

f 

«n 

ft  <  * 

«o 

cm  a  x  z 

a 

* 

W 

>  to 

z 

in 

X 

<  ►  o 

H 

t— 

w 

> 

-  z>  a. 

1— 

< 

< 

“>  x  > 

z> 

z  a 

•t 

t~  a  < 

a 

o  o 

f 

x  to  * 

z 

LO 

'O 

in 

>  >  X 

X 

z 

f 

C\J 

<  *  < 

o  o 

X 

•  r>  a 

UJ  u 

ft 

> 

3  X  > 

z 

1- 

z 

•  a  < 

ft 

< 

o 

X  CO  • 

1- 

3: 

D 

ft 

m 

•  X 

3 

Z 

cO 

W 

Z  D  U 

o 

• 

t— 

z 

<  o  x  a 

a 

-1 

1- 

a 

Q  x  a  x 

CE 

z 

X 

X  X  •>  > 

O 

• 

o 

*—* 

<  o  <  < 

CO  3 

u 

o 

a> 

rH 

ft  CM 

■)  C  IL  IL  D 

N  N  D  N  J 

vj  M  tf-  -  ■  ■*■ 


i 


i 


I 


CONTINUED 


CONTINUED 


140  SUMDF  =  SUMDF+DF 


ooO'CHrvjri'tir.  iCNeoO'  o<Hrg«ci4miOf~-ooo>o<-if\j<n.$ 

4,4,mtnininmir>if\inmini0i0>0'0'0<0'0v0'0'0f,“f^t,~f^f** 

COCOOOCCCCCOCOCCCOOOCCOC  COC 


•  c 

a  3 

m 

s  X 

CM 

<  UJ 

w 

« 

m  > 

X 

1- 

u  ► 

V 

< 

» 

*  c 

1- 

s: 

0 1- 

z 

"> 

• 

3  * 

Uj 

* 

DC 

X  3 

UJ 

< 

U  < 

O' 

in 

K 

1- 

U 

V 

<NJ 

in 

K  •* 

a 

1- 

w 

1- 

C.  x 

\ 

z 

D 

» 

3  X  I 

X 

h- 

UJ 

2. 

oc 

X  X 

CO 

Z 

U 

Lu 

< 

UI  3 

« 

UJ 

CL 

cO 

K 

»  in 

(M 

u 

U 

> 

in 

X  • 

• 

CL 

Q 

• 

uj 

U  fM 

»o 

UJ 

» 

uj  in 

LL 

a 

X 

m 

UJ 

x  * 

•> 

00 

CM 

5: 

•> 

L;  1  rH 

CL 

i 

•> 

w 

*— 1 

a 

*  X 

O 

oc 

CO  — 

Zj 

h— 

> 

U-  < 

O' 

• 

•  in 

2: 

» 

< 

0  X 

CL 

(M 

cn  — 

nj 

LU 

3 

« 

»  3 

Uj 

• 

rH  • 

1— 

to 

» 

5 

<  * 

00 

u  o 

3 

*  * 

Z 

>  *-  x 

_J 

u 

•  UJ 

a 

U' 

3  < 

< 

• 

oc  ;-j 

z 

m  — 

► 

•> 

uj  x 

■s  a:  — 
(DO® 

•—  O'  •— 

►  oc  • 

5.  uj  5 

o  o 
a  o  a 

Ul  z  U 


u  a  u. 
x 

u  <  u 
O  to  o 


in  C 
c  *  x 
o  o  c 

•  vO  <N 
C  •— 1  — • 
II  I- 

<  X  < 

* -  u  x 

J  2  K 
IUDO 
O  Q.  U. 


K  ? 

*-  O  • 

3  a  _j 

Q  a  cd 

Z  Z  3  • 

UJ  *->  lO  2 

t- 

cc  a  uj  a  in 

O  3  3  C  «■* 

o  uj  u  a 

NNONJ 
fsj  M  *  K  *  U 


uj  *  m  a  >  c.  3 

Z  m  nj  <  »  ►  ► 

»-«  r\j  ■—  I—  <  Q.  x 

I  - Cl  \T,  C,  X  < 

3D>ZJ<i 

o  x  <  uj  <  m  \l 

a  uj  *  »  *>  • 

•  —  ->  X  <  -u 

i—  in  •—  c  •  * 

DU)Ma  i  zh 

a  •  —  >  ->  u,  x 

Z  in  x  <  -  O  < 

—  ~  >  *  3  a.  x 

<  <  ”>  x  >  “> 

I I  •  »  •— i  Li_  <  • 

o  -*  x  m  •  a. 

«o  m  >  >  x  X 

Z  —  CM  <  •  <  < 

c  x  —  •  3  cl  in 

—  >  26  X  >  - 

K  Z  *  •  U-  <  * 

<  o  —  x  m  »  o: 

3  •-  m  •x*- 

z  m  ~  z  3  u  * 

—  z  <  o  rr  uj  x 

h  U  O  I  UJ  s  < 

Z  X  X  X  -  >  X 

O  •—  <  O  <  < 

u  a  u 


CONTINUED 
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Figure  4.  Continued. 
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Figure  4.  Continued. 
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Figure  4.  Continued. 
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Figure  4.  Continued. 
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Figure  4.  Continued. 
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Figure  4.  Continued. 


o  h  n  n  ^  iftior^ooa  OHN(n4if'»ONooO'OHrg(fij'if\ 
O'  O'  O'  O  O  O'  O'  O  O  O'  OOGGCOOOOOr— |>— ll— irHr-irH 

>?«?«*  <1  ^  ^  4  4,<f>4'inininininininininininininininin 
coooocoocooccccoccooccoooc 


m 

x 

< 

x 


in 

rvj 


u 

•> 

ce 

< 

t- 

LQ 


>H  • 

•  C 

a  x 

X  X 
<  UJ 
m  > 

UJ  • 

•  c 

C  I- 

3  • 

2:  x 

Ll  < 

•  (- 
t-  • 


1/5 

X 

u 

lu 

X 

UJ 

• 

in 


X 

* 

X  X 

a 

no 

x 

or 

X  X 

X 

iD 

u: 

< 

Uj  X 

< 

0> 

m 

t— 

•  tn 

in 

> 

m 

x  • 

lL 

I— 

• 

IX 

U  nj 

w 

— 

#• 

Uj  m 

Ul 

X 

m 

u. 

X  • 

c 

> 

X 

•Ll  h 

< 

w 

•—4 

CL  •  X 

X 

X 

K 

>  Ll  < 

m 

^1. 

X 

m 

<  O  X 

•— • 

rH 

UJ 

X 

•  •  X 

0. 

w 

in 

m 

• 

X  <  • 

X 

a 

o 

UJ 

• 

z 

>  h-  X 

< 

> 

z 

z 

UJ 

<  X  < 

in 

X 

< 

•— < 

m 

m 

•  Uj  X 

Ll 

< 

•» 

in 

►— 

• 

m 

QC 

>  a  _j 

c, 

X 

— . 

X 

m 

CsJ 

< 

»  •  » 

•> 

-J 

-J 

•— > 

o 

(\l 

t— 

<  a  x 

X 

* 

» 

w 

or 

a 

m 

o  x  < 

a 

rH 

CD 

D 

> 

z 

X  <  2 

X 

II 

II 

— 

5  — 

X 

5: 

< 

UJ 

<  in  ic 

< 

rH 

00 

CO 

UJ  GO 

i/J 

LL 

* 

•  >  * 

in  c  . 

_l 

• 

•» 

• 

<r  • 

r 

X 

X  <  X 

>  o 

* 

>  ^ 

m 

— ‘ 

in 

*— « 

*  #• 

*«— 

rH 

rH 

+  *H 

o 

vC 

<\ 

a 

X  z  — 

I  c 

rH 

UJ 

Uj 

-*  LU 

z 

• 

W 

> 

x  m  x 

n 

u 

-J 

in 

►H 

in 

•  in 

< 

in 

A*. 

< 

•  o  < 

X 

UJ  — 

•* 

* 

» 

w  #• 

| 

rH 

in 

w 

> 

x  a  x 

< 

r>  rH 

-J 

in 

X 

in 

D  m 

in 

< 

< 

X 

X  >  X 

m  Ll 

w 

M 

X 

*— 4 

ZL  -* 

! 

c 

u. 

e> 

* 

M 

Lu  <  * 

>  a  x 

< 

U! 

#1 

UJ 

. 

z 

O 

— * 

X 

m  •>  a 

* 

•  u 

X 

w 

m 

w 

m 

</)  m 

< 

in 

> 

>  X  X 

■a- 

x  t- 

< 

•> 

*— • 

» 

•>  »—i 

1 

in 

z 

z 

•w 

no 

< 

*  <  < 

X 

X  — 

X 

* 

V 

•* 

^  * 

' 

c 

O  X 

«— 

•> 

x  a  o', 

U 

U  3: 

-J 

r> 

r 

m 

♦  v0 

UJ 

. — 1 

in 

> 

:< 

X  >  — 

u 

uj  m 

+ 

►H 

►— * 

i 

z 

t- 

X 

z 

» 

U-  «  • 

X 

z~ 

rH 

•i 

X 

*> 

I 

•  X 

►- 

< 

o  o 

X 

in  *  or 

> 

>  LU 

II 

c 

fVJ 

c: 

nj 

ro 

1- 

X 

u 

r— < 

m 

•  x  — 

• 

m 

_! 

'O 

— 

fH 

w 

OC  w 

! 

H 

H 

X 

z 

m 

W 

z 

x  u  • 

X 

z  z 

1— 

h- 

z 

i 

m 

m 

o 

t — • 

*H 

z 

< 

o 

X  UJ  x  u 

O  Lu 

O 

X 

< 

X 

< 

X  < 

QC 

! 

a 

Ct 

oc 

t- 

3:  UJ 

o 

X 

lu  x  <  uj  x  m 

in 

u 

s: 

X 

^  r 

X 

! 

z 

z 

CD 

z 

ZL 

5: 

X 

•  >  X 

a: 

z 

(T 

2  a 

1- 

o 

< 

< 

X 

o 

• 

•— « 

< 

o 

<  <  — 

>  O  U. 

O 

D 

O 

X 

O 

3  O 

U! 

; z ■ 

m 

m 

m 

u 

> 

O 

u 

U  I-H 

o  a 

Ll 

a 

U 

a  a. 

QC 

UJ 

*- 

X 

1 

2 

3 

<3- 

m  a  uj 

QC 

in 

o 

o 

o 

o 

o 

O 

O  X  _J 

O 

*— » 

m 

sO 

00 

O' 

->  o  uj 

Ll 

C5 

N  N  C5 

isl 

_l 

rsj  M  * 

NJ 

IX 

66 


CONTINUED 


Figure  4.  Continued. 
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Figure  4.  Continued. 
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Figure  4.  Continued. 
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Figure  4.  Continued. 
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Figure  4.  Continued. 
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Figure  4.  Continued. 
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Figure  4.  Continued. 
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GLOSSARY  OF  SYMBOLS 


Cov  (a,  b) 


5 

Sum  of  all  the  L  elements  of  the  last  inverse 
matrix  for  sample  i  in  the  iterated  least-squares 
procedure 

Expected  value  of  the  background  for  the  jtfl 
observation  on  the  i'^  sample. 

Counts  observed  during  period  T,  Tjj. 

Counts  of  background  activity  observed  during 
period  T*,  Tjj. 

Estimated  covariance  of  a  and  b. 

Number  of  degrees  of  freedom  associated  with 
the  historic  estimate  of  the  sampling  error,  aso. 

The  expected  value  or  arithmetic  average  value 
operator,  such  that  in  the  discrete  case 

1  " 

E(y)  =  lim  -  V  yh  ,  where  y  is  any 
n-°°  "  h-1 

random  variable  and  y1(  y2,  V3,  are  values 
that  y  takes  on  in  successive  occurrences. 

Number  of  degrees  of  freedom  associated  with 
the  historic  estimate  of  the  mechanical  error, 


Number  of  degrees  of  freedom  associated  with 
the  current  estimate  of  the  mechanical  error, 
am'  ,  based  on  the  historic  estimate  amo . 


ft  =  £  (Ji  -  L) 


Student's  t  probability  density  function. 

An  integer,  h  =  0,  1,  2,  ... 

Number  of  aliquot  samples  of  radioactive 
solution. 

Index  to  label  sample;  i  =  1,  2,  3, ..  ,  I. 


N*>  Njj 


Number  of  observations  of  the  i**’  sample. 

Index  to  label  observation;  j  =  1,  2,  3 . Jj. 

A  normalising  constant  for  the  i’*1  sample!  on  the 
k  th  iteration  in  the  computation  of  the  weights 
Wjjk  so  that 
Ji 

£  wijk  =  1 

i  “  1 


Iteration  index  used  in  least-squares  curve¬ 
fitting  procedure;  also  in  estimating  os 

Number  of  radioactive  isotopes  in  the  source 
solution. 

Index  to  label  isotope;  C  =  1,  2,  3, ....  L. 

The  mean  of  a  normally  distributed  random 
variable.  Used  for  illustrative  purposes  only 
when  discussing  the  mechanical  error. 

Counts  per  minute  from  both  sample  and  back 
ground  corrected  for  deadtime  of  the  detector. 


Apparent  counts  per  minute  from  background 
activity. 

C*  C*: 

N*  =  —  and  Njj  = 

ij 

Probability. 

The  mean  of  a  normally  distributed  random 
variable.  Used  for  illustrative  purposes  only 
when  discussing  the  sampling  error. 

A  normally  distributed  random  variable  with 
mean  zero  ana  standard  deviation  om  labeled 
for  the  j,fl  observation  of  the  i,h  sample  used 
in  the  simulation. 


1 


I 


Number  of  observations  of  the  i**'  sample. 
Index  to  label  observation,  j  =  1,  2,  3, Jj. 


A  normally  distributed  random  variable  with 
mean  zero  and  standard  deviation  os  labeled 
for  the  i^'  sample  used  in  the  simulation. 


A  normalizing  constant  for  the  i**1  sample  on  the 
kth  iteration  in  the  computation  of  the  weights 
Wjjk  so  that 
Jj 

Z)  wijk  =  1 
i  =- 1 

Iteration  index  used  in  least-squares  curve¬ 
fitting  procedure;  also  in  estimating  as 

Number  of  radioactive  isotopes  in  the  source 
solution. 

Index  to  label  isotope,  C  =  1,  2,  3 . L. 


For  the  k,h  iteration  the  weighted  sum  of 
squares  of  deviations  of  the  observed  counts 
per  minute  for  the  i'*1  sample  from  the  cor¬ 
responding  fitted  curve  divided  by  the  degrees 
of  freedom 

1  Ji  ~  2 

®ik  ’  j.  _  "l  Z  wijk  "  F*.ik  ^ij> 


Estimated  variance  of  p  ■.  . 


Time  in  minutes  expended  in  observing  a  sample 
of  radioactive  solution. 


The  mean  of  a  normally  distributed  random 
variable.  Used  for  illustrative  purposes  only 
when  discussing  the  mechanical  error. 

Counts  per  minute  from  both  sample  and  back¬ 
ground  corrected  for  deadtime  of  the  detector. 


T*  T* 

1  >  1  ij 


■•fa 


and  N 


Apparent  counts  per  minute  from  background 
activity. 

C*  cii 

N*  =  — -  and  N-  =  — — 


Probability. 

The  mean  of  a  normally  distributed  random 
variable.  Used  for  illustrative  purposes  only 
when  discussing  the  sampling  error. 

A  normally  distributed  random  variable  with 
mean  zero  and  standard  deviation  om  labeled 
for  the  jth  observation  of  the  i,h  sample  used 
in  the  simulation. 


Var  (a) 


Time  in  minutes  expended  in  observing  the 
background. 

Age  of  sample  in  minutes  ai  the  beginning  of  the 
observation  period  T  or  T Sj .  where  t  and  tjj  are 
reckoned  from  a  historic  instant  of  time  when 


A  historic  instant  of  time  when  the  age  of  all 
samples  used  in  a  particular  experiment  are 
reckoned  to  be  zero  minutes  old;  that  is,  t  =  0. 

Variable  defined  by  Student's  t  density  function. 

If  a  is  the  number  of  samples  for  which  a  par¬ 
ticular  background  reading  was  ued  in  the 
compc  ations,  and  b  is  the  number  of  times 
it  was  used  just  for  the  i,h  sample,  then  Uif  =  ab. 

Weight  for  the  i’*1  sample  on  the  k1*1  iteration 
used  to  determine  p  0  and  os . 

Estimated  variance  of  a. 

Weight  for  the  j1*1  observation  of  the  i’^1  sample 
on  the  i  ,h  iteration. 
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GLOSSARY  OF  SYMBOLS  (Continued) 


X(t,T;Xg) 


Independent  variable  in  the  iterated  linear  least- 
squares  curve  fit  (corrected  for  decay). 


X(t,T;Xg)  =  e 


xf 1  -  e  ® 


XgT 


X 


Ei) 


Same  as  X(tjj,  T jj ;  Xg ). 


Individual  point  estimate  of  the  variation  due 
to  the  Poisson  phenomenon  for  the  j1^1  obser¬ 
vation  of  the  ith  sample. 


Crude  estimate  of 


fin 

0. 

T. 


Individual  point  estimate  of  the  variation  due 
to  mechanical  causes  for  the  jth  observation 
of  the  ith  sample. 


Crude  estimate  of 


Crude  estimate  of 


Random  error  on  the  j,h  observation  of  the 
i1*1  sample. 


4.{t) 


p.  (t> 


/ig«(t.T) 


Expected  instantaneous  counts  per  minute  from 
isotope  L  of  sample  i  at  age  t,  corrected  for 
background  and  deadtime  of  the  detector 

'11  °  * 

pgilt)  =  PCj  e 


Expected  value  ol  pg.  (t);  that  is, 

1  1 

fig  (t>  =  lim  j  ^2  Pgj(t) 


I  -»°° 


i  =  1 


Same  as/ig  (t)  except  for  all  isotopes: 

L 

/■c*1* =  53  <*> 

£=  i 

Expected  counts  per  minute  for  a  period  T  from 
isotope  8  of  sample  i  at  age  t,  corrected  for 
background  and  deadtime  of  the  detector. 


M.j  (t,  T) 


Same  aspgf  (t,  T)  except  for  all  isotopes  in  the 
sample;  that  is, 


M.i(t,T)  =  £  MgjttiT) 


T/,T?ij 


«k 


h 


Assumed  counts  per  minute  arising  solely  from 
the  sample.  It  is  also  the  dependent  variable 
in  the  regression  analysis. 

tj  =  N  -  N*  and  rjjj  =  Njj  -  Njj 


Pg  (t.  T) 


The  population  mean  of  pg.  (t,  T);  that  is, 

pg  (t,  T)  =  lim  -j-  53  A>gj(L  T) 
l-00  i«1 


A  normalizing  constant  in  the  estimation  of 
the  pg  's  such  that 

53  vik  = 1 

i= i 


P..(t,T) 


The  population  mean  of  p  j  (t,  T);  that  is, 

p  (t,  T)  =  l|m  T  53  ^  i<l'T* 

I  -  i  =  i 


A  decay  constant  in  reciprocal  minutes  for  the 
8th  isotope,  fi=  1,  2, ....  L.  It  is  assumed  known 
and  free  from  error  in  this  report. 


■  E  “eA" 
s- 1 


A 


4<  t) 


4W 


H[ .  It) 


Mg,  <t,  T) 


Expected  instantaneous  counts  per  minute  from 
isotope  L  of  sample  i  at  age  t,  corrected  for 
background  and  deadtime  of  the  detector: 

'  /.i  0  "^£l 

pgi(t)  =  Pg,  e 

Expected  value  of  pg.  (t);  that  is, 

1  1 

jig  (t)  =  lim  y  £  pg,(t) 
l*°°  i-i 

Same  aspg  It)  except  for  all  isotopes: 

L 

Ji*.<t)  =  X  **£  It) 

6-1 

Expected  counts  per  minute  for  a  period  T  from 
isotope  £  of  sample  i  at  age  t,  corrected  for 
background  and  deadtime  of  the  detector. 


M,i  (t.  T) 


Pg  It.  T) 


P..(t,T) 


Same  as  pg.  |t,  T)  except  for  all  isotopes  in  the 
sample;  that  is, 

L 

M.ilt.T)  =  E 
£■1 

The  population  mean  of  pg,  It,  T);  that  is, 

pg  It,  T)  =  lim  ^  Pgi  It*  D 


The  population  mean  of  p ,  (t,  T);  that  is, 

p  (t,T)  =  lim  J  X  M.jlt.T) 


|  -♦  OO 


i  =  1 


=  x  Mg  it.  T) 

£  =  t 


^£i 


J'.i 


O 


^£. 


o 


M. . 


Crudely  the  expected  counts  per  minute  at 
time  t  =  0  from  isotope  £  of  sample  i: 

Pg!  =  lim  Pgi  (0,  T) 

T  -*  OO 


Same  as  pg.  except  for  all  isotopes  of  sample  i: 

L 

O  r*^  O 

^  =  E  ^£i 

£'i 


The  population  mean  of  /-tg . ;  that  is, 

0.1  O 

=  "m  J  E  ^£i 


I  -  OO 


i- 1 


Same  as  pg  except  for  all  isotopes: 


O  o 

^ =  E  n. 

£-t 


For  the  jth  observation  of  the  ith  sample  of  the 
generated  data,  an  intermediate  quantity  (counts 
per  minute)  that  includes  the  effects  of  both 
sampling  and  mechanical  errors  but  is  free  of 
Poisson  effects. 

$ij  =  0  +  ri><1  +  Qjj)  Ji.i  (tjj,  Tjj) 

The  Poisson  parameter  T,j(£,j  +  B,|)  is  required 
to  generate  C,j. 

Used  in  general  to  indicate  a  variance  (may  be 
porportional  in  some  cases).  For  example  oa2 
is  used  to  designate  the  variance  of  the  quantity 
a.  Likewise,  ca2  would  be  an  estimate  of  the 
variance  of  a. 


Proportional  mechanical  error;  that  is,  a 
proportional  error  due  to  all  causes  other  than 
sampling  error  and  errors  that  result  from  the 
Poisson  phenomenon. 
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OF  8YMB0LS  (Continued) 


An  estimate  of  am  on  the  k^1  iteration  of  the 
least-squares  process  is  a  historical  esti¬ 
mate  required  for  startinq  the  iterated  least- 
squares  process. 


The  appropriate  historic  estimate  of  the 
proportional  mechanical  error  based  on  f0 
degrees  of  freedom. 

The  current  proportional  mechanical  error 
based  on  o^. 

An  estimate  of  am-  on  the  kth  iteration. 

Proportional  sampling  error. 

An  estimate  of  at  on  the  kth  iteration  of  the 
least-squares  process. 

The  appropriate  historic  estimate  of  the 
proportional  sampling  error  based  on  d 
degrees  of  freedom. 


The  current  proportional  sampling  error  based 
on  o*, . 


An  estimate  of  of<  on  the  kth  iteration. 


Total  estimated  variance  of  p .  , 


2 


°~o 

H.. 


Estimated  variance  of  p°  due  to  mechanical 
effects. 


Estimated  variance  of  p  due  to  Poisson  causes. 


Estimated  variance  of  p  due  to  sampling 
effects. 

Deadtime  of  the  detector  in  minutes. 


o  When  used  as  a  superscript,  it  indicates  the 

quantity  essentially  at  time  t  ■  0;  for  example, 


When  used  as  a  subscript,  it  acts  as  a  placeholder 
for  that  particular  subscript.  Interpretation  of 
usage  is  usually  different  in  each  case. 

A  bar  placed  over  any  symbol  indicates  an 
arithmetic  average  or  the  mean  of  a  finite 
number  of  items;  for  example, 

1  n 

n  £  * 
h-1 

A  tilda  placed  over  any  quantity  indicates  it  is 
an  estimate  of  the  true  quantity;  for  example, 
o 

Pj(  is  an  estimate  of  pg, . 
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